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ABSTRACT 

The nature of the interaction between galaxies and the intergalactic medium (IGM) is one of the most funda- 
mental problems in astrophysics. The accretion of gas onto galaxies provides fuel for star formation, while 
galactic winds transform the nearby IGM in a number of ways. One exciting technique to study this gas is 
through the imaging of hydrogen Lya emission. We use cosmological simulations to study the Lya signals 
expected from the growth of cosmic structure from z — 0-5. We show that if dust absorption is negligible, 
recombinations following the absorption of stellar ionizing photons dominate the total Lya photon produc- 
tion rate. However, galaxies are also surrounded by "Lya coronae" of diffuse IGM gas. These coronae are 
composed of a combination of accreting gas and material ejected from the central galaxy by winds. The Lya 
emission from this phase is powered by a combination of gravitational processes and the photoionizing back- 
ground. While the former dominates at z 0, collisional excitation following photo-heating may well dominate 
the total emission at higher redshifts. The central regions of these systems are dense enough to shield them- 
selves from the metagalactic ionizing background; unfortunately, in this regime our simulations are no longer 
reliable. We therefore consider several scenarios for the emission from the central cores, including one in which 
self-shielded gas does not emit at all. We show that the combination of star formation and cooling IGM gas can 
explain most of the observed "Lya blobs" at z ^ 3, with the important exception of the largest sources. On the 
other hand, except under the most optimistic assumptions, cooling IGM gas cannot explain the observations on 
its own. 

Subject headings: cosmology: theory - galaxies: formation - intergalactic medium - diffuse radiation 



L INTRODUCTION 

The distribution of matter on large scales is one of the key 
ingredients of any cosmological paradigm. A number of ob- 
servations constraining the overall cosmological paradigm, 
together with numerical simulations, have shown that the 
"cosmic web" provides a good description of the Universe on 
the largest scales. In this model, structure formation proceeds 
through the collapse of matter into sheets, filaments, and even- 
tually galaxies and galaxy clusters. However, on small scales, 
baryons and dark matter are subject to different forces, and 
the resulting complications have severely limited the success 
of ab initio models in predicting the structure and evolution 
of galaxies. The interactions between galaxies and the inter- 
galactic medium (IGM) are particularly fundamental to our 
understanding of the formation and growth of galaxies. The 
accretion of gas onto galaxies provides fuel for star forma- 
tion, while feedback processes (both radiative and mechani- 
cal) from star formation can transform the IGM. 

Quasar absorption lines are the most well-known way to 
study the IGM at moderate redshifts, and there is no doubt 
that the Lya forest has provided a wealth of insight into the 
IGM, the distribution of gas around galaxies, and the nature 
of the interactions between the two (see Rauch 1998 for a re- 
view). An important aspect of the forest is that absorbers are 
localized in redshift so that each can be easily separated in the 
redshift direction. However, because the forest probes only 
isolated lines of sight, it has proved difficult to determine how 
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the absorption systems relate to the cosmic web and to galax- 
ies. These are particularly intriguing questions because the 
answers will shed light on the formation and growth of galax- 
ies. An alternate approach is to search for Lya emission from 
hydrogen in galaxies and in the IGM. This method has the 
advantage of permitting a direct reconstruction of the three- 
dimensional gas distribution, which will allow much cleaner 
constraints on the interplay with galaxies. Moreover, the Lya 
radiation directly probes an important fraction of the cool- 
ing radiation at the temperatures relevant for the formation 
of galaxies, and it is often more readily detected than contin- 
uum radiation. For these reasons, this approach has received 
a large amount of attention, despite the challenging nature of 
the observations. 

In particular, narrowband Lya-selected surveys have 
emerged in recent years as an efficient way to identify 
and localize high-redshift objects. Such surveys have been 
used to detect faint galaxy associations at moderate red- 
shifts (Steidel et al. 2000; hereafter SOO), high redshift 
galax ies (iHuet al.ll2002t IKodaira et ab 2003; R hoads et aTl 
[2004^, the filamentary structure of galaxies at high redshifts 
(Moller & Fvnbo 20011), th e host galaxies of Lyma n-limit and 
damped Lya systems (e.g., Fvnbo et al. 199?, "2000), and dif- 
fuse emission near radio galaxies (McCarthv et al. 1987). On- 
going efforts to detect Lya emission from high-column den- 
sity absorbers in the IGM ha ve already yiel ded interesting 
constraints ( Francis & Bland- Hawthornll2004ft and will soon 
reach the predictions of simple models (IGould & Weinberg I 
1996). 

The last effort is an example of one of the most exciting 
possibilities for Lya surveys: imaging gas outside of galax- 
ies. To date, most Lya-selected sources are compact star- 
forming galaxies. Ionizing photons from young, massive stars 
are absorbed by the interstellar media of the galaxies and re- 
radiated (following recombinations) as Lya photons. How- 
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ever, the ubiquity of hydrogen also permits the direct de- 
tection of the IGM. One example is the diffuse Lyg emis - 
sion often found near radio galaxies (^ McCarthy etaTHToSTl) . 
thought to be shock-heated gas interacting with the radio jet. 
Another, perhaps more representative, set of objects are ex- 
tended "Lya blobs" at z ^ 3 with little or no continuum 
emission. SOO detected two distinct regions of diffuse Lya 
emission around (but not centered on) three of the Lyman- 
break galaxies in their field (which contained 27 such galax- 
ies). These objects have Llyq ^ 10'*'* erg s^* and physical ex- 
tents of ~ kpc. More recent observations revealed 33 
smaller blobs in an overlapping field, many clustered around 
known galaxies ('Matsuda et al. 2004'; hereafter M04). These 
Lya blobs are particularly intriguing because they sit at the 
interface between galaxies and the IGM and should probe the 
interactions between the two. 

There are a number of processes that can cause Lya 
emission from the IGM, most of which involve interactions 
with galaxies. Hogan & Wevmann ( 1987) suggested that 
optically thick, high-column density systems should absorb 
ionizing photons from the metagalactic radiation field and 
re-radiate a substantial f raction of them as Lya photons. 
iGould & WeiTlto ? (TOQ^ refined this argument and showed 
that deep observations with 8-meter class telescopes could de- 
tect the resulting emission. Interestingly, the luminosity of 
a system scales with the local ionizing background in this 
scenario (so long as it remains optically thick). Thus the 
Lya blobs could be extreme examples of this effect, in which 
ionizing photons from luminous, embedded sources are re- 
processed into Lya photons. The largest blob does con- 
tain a luminous submm (though optically invisible) source 
with an inferred star f ormation rate (SFR) > 500/;^^ Mq yr^' 
iChapman et alJ2()(H ). which could power the extended emis- 
sion. An alternate possibility is the assembly of gravita- 
tionally bound objects. In the traditional view, baryons are 
shocked as they accrete onto halos (White &Rees 1978). 
They then shed their gravitational energy through radiative 
cooling, and a large fraction of the energy could be lost 
as Ly a radiation (e.g., Haiman et al. 2000). Fardal et ^ 
(|2001) examined this mechanism with cosmological simula- 
tions and argued that it could explain the luminous Lya blobs, 
although they find weaker shocks than e xpected (see also 
iBirnboim & Dekel''2003^ 'Keres et al."2004V Finally, there is 
eviden ce fo r strong galactic winds near the m ost luminous 
blobs ('Steidel et al. 2000; Taniauchi et al. 20()T; "Bow er etalJ 
1^04), and the Lya emission could be a result of the interac- 
tion of this wind with the surrounding IGM. Similar interac- 
tions probably explain the diffuse Lya emission surrounding 
radio galaxies at high redshifts (McCarthy et al. 1987). 

Unfortunately, these scenarios are difficult to distinguish, 
and there has been relatively little theoretical investigation of 
the possibilities. In Furlanetto et al. ( 2003, hereafter F03), we 
used a state-of-the-art cosmological simulation to predict the 
Lya emission from IGM gas at z < 0.5. We found that most 
galaxies were surrounded by "coronae" of Lya emission a few 
times larger than the embedded galaxies. We argued that the 
emitting gas is slowly cooling onto the galaxies, so that the 
Lya emission traces the growth of galaxies. In this paper, we 
extend our work to higher redshifts < z < 5 and examine in 
detail the origin of the emission and of the gas responsible for 
it. We show that Lya emission can be particularly powerful 
during the era at which galaxy assembly peaks (z ^ 1-3) and 
we discuss the relative importance of the different physical 



mechanisms that are responsible for the production of Lya 
photons from diffuse gas. 

At the same time, it is useful to consider Lya emission from 
inside of galaxies. Massive, hot stars produce ionizing pho- 
tons, most of which are absorbed by the surrounding hydro- 
gen gas. As this gas recombines, it radiates Lya photons; the 
Lya luminosity thus offers a (rough) measure of the star for- 
mation rate, although dust and radiative transfer through the 
galaxy and IGM can dec ouple the two in individual objects 
(e.g., Shapley et al. 2003). The SFR in turn depends (at least 
crudely) on the rate at which gas accretes onto the galaxy, so 
this mode also contains information on how galaxies grow. 
The stellar ionizing photons can also escape into the IGM, il- 
luminating the neighborhood of each galaxy and boosting the 
brightness of the Lya coronae. We thus also use the simu- 
lations to predict the expected emission that can be traced to 
stellar ionizing photons. We will show that the emission from 
the IGM and from star formation both contribute to any given 
source, but that the former is more spatially extended. 

We first discuss the numerical simulations and our analysis 
procedure in |2] We then describe how to compute the Lya 
emissivity of a parcel of gas in |3l We review the properties of 
the IGM, and how they relate to the Lya emission, in |3 We 
then describe our results in terms of the statistical properties 
of maps in ^and in terms of individual objects in ^ Finally, 
we conclude in ^ 

We present most of our results in terms of surface bright- 
ness. We use two (interchangeable) sets of units; for refer- 
ence, 1 photon s~' cm^^ sr^' ~ 9.6 x 10~^-'[4/(l +z)] erg 
s^' cm^^ arcsec^^. The redshift factor appears because 
Aobs — Aq(1 +z), where Aq = 1215.67 A is the rest wave- 
length of the hydrogen Lya transition. 

2. SIMULATIONS AND ANALYSIS PROCEDURE 

We use the suite of cosmological simulations described 
by Springel & Hernquist ( 2003b) and Nagamine et al. ( 2004). 
These smoothed-particle hydrodynamics (SPH) simulations 
include a multiphase description of star formation incorpo- 
rating a prescription for galactic winds ( Springel & HernauisJ 
'2003a) and were performed using a fully conservative form 
of SPH (Springel & Hernauist 2002). They also accounted 
for the presence of a uniform ionizing background as de- 
scribed by Haardt & Madau (1996). The simulations as- 
sume a ACDM cosmology with fi,,, = 0.3, f^A = 0.7, fit, — 
0.04, //,) = 100/! km s-i Mpc"' (with h = 0.7), and a scale- 
invariant primordial power spectrum with index n ~ I nor- 
malized to (jg ~ 0.9 at the present day. These parameters 
are c onsistent with the rn ost recent cosmological observations 
(e.g.,|SEerieretal]200l. Because we wish to study a variety 
of redshifts and physical scales, we will use several different 
simulations. We summarize their main parameters in Table|2| 

To make maps of the Lya emission, we wish to compute the 
surface brightnes s of a slice of the simulati on with fixed z, Az, 
and angular size. iFurlanetto et alJ ( l2004aft describe our anal- 
ysis procedure in detail, so we only summarize it here. We 
first randomly choose the volume corresponding to this slice 
from the appropriate simulation output. We then divide the 
projected volume into a grid of A^^ pixels; in most of what fol- 
lows we choose = 512. For each particle in the simulation, 
we check whether it lies within the slice of interest. If so, we 
compute its Lya emissivity by linearly interpolating (in log 
space) the grids described in ^ assuming that the particle 
is cubical with a uniform density (this is reasonable because 
luminous particles are much more compact than our map pix- 
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Table 1 . 



Name 


L 


Resolution 


"Jgas 


e 




Winds 


Q5 


10.00 


2 X 324^ 


3.26 X 10' 


1.23 


2.15 


Strong 


Q4 


10.00 


2x216^ 


1.10 X 10* 


1.85 


2.75 


Strong 


Q3 


10.00 


2 X 1443 


3.72 X 10* 


2.78 


2.75 


Strong 


P3 


10.00 


2 X 1443 


3.72 X 10* 


2.78 


2.75 


Weak 


03 


10.00 


2 X 1443 


3.72 X 10* 


2.78 


2.75 


None 


D5 


33.75 


2 X 3243 


1.26 X 10' 


4.17 


1 


Strong 


G6 


100.0 


2 X 486^ 


9.67 X 10' 


5.00 





Strong 


G5 


100.0 


2 X 3243 


3.26 X 10** 


8.00 





Strong 



Note. — All simulations are described in Sorinael & Hernauist 1 2003lJ, 
with the exception of G6 (which is described in Naeamine et al. 2004). Here 
L is the box size in /j^' comoving Mpc and e is the softening length in /i^' 
comoving kpc. The "Resolution" column gives the initial number of particles 
(including l^oth gas and dark matter). The particle mass mga, is listed in 



h ' Mq. "Strong" and "weak" 
242 km s~', respectively. 



winds have nominal velocities of 484 and 



Thick: photoionized 
Thin: Lyy pumping 
Dot-dashed: CTE 




Fig. 1. — Lya emissivities at z = 3. The thick curves show ea/n^ for 
log(«H/cm^^^) = —2, —4, and —6 (solid, dashed, and dotted curves, respec- 
tively). The thin curves show the pumping emissivity from Ly7, assuming an 
optically thin medium, for the same densities. The dot-dashed curve shows 
the Lya emissivity for gas in collisional ionization equilibrium (CIE). 



els). We then determine which pixel(s) the particle overlaps 
and add the particle's surface brightness to them, weighting 
by the fraction of the pixel subtended by the particle. Finally, 
we smooth the maps using a Gaussian filter with a width of 
four pixels; the resolution /S.6 that we quote is the FWHM of 
this smoothing function. 

3. LYcj EMISSION 

Our emissivity calculations follow, for the most par t, the 
same procedures as in FOB and lFurlanetto et al.l i2004aft . We 
first divide the gas into three regimes: optically thin, self- 
shielded, and star-forming. The first component consists of 
gas elements that are optically thin to the background ioniz- 
ing radiation field. We u se the photoionizing backgrounds of 
iHaardt & Madaul ( l200l . which include emission from galax- 
ies and quasars as well as reprocessing by the IGM. We de- 
note the total ionizing rate by r = ri2 x 10^'^s^^ For ref- 
erence, we have Tp = (0.086,0.790, f.50, 1.15,0.761,0.565) 
atz= (0,1,2,3,4,5). 

3.1. Optically Thin Gas 

To compute the Lya emissivity of optically thin gas, we 
construct a grid in density ne with spacing Alogne — 0.25 
in the range —6 < log«H < — 1 and Alog/ie = 1.0 elsewhere 
(here hh is measured in cm^^) and another grid in temper- 
ature with AlogT = 0.05 in the range 3.5 < logT < 5 and 
A log T ~ 0.25 elsewhere (here T is measured in degrees K). 
We then use the Cloudy 96 photoionization code (beta 5; 
lFerlandl[2003l) to compute the emissivity at each of the grid 
points, assuming solar metallicity. (The results are insensitive 
to the metallicity, so long as it is near or below this value.) 

The thick curves in Figure [2 show Eq/wh fo'" three differ- 
ent densities at z = 3 in the optically thin case; other redshifts 
are qualitatively similar. At small riu, collisional excitation 
can be ignored and Ea essentially follows the temperature de- 
pendence of the recombination coefficient. In this regime, 
the emissivity is independent of the amplitude of the ioniz- 
ing background. As the density increases, collisional pro- 
cesses become increasingly important and cause the increase 
in ea/nji alT 10'*^-10^ K. We note that in this regime exci- 
tation dominates the emissivity, with recombinations (follow- 
ing either collisional or photoionization) making only a small 
contribution. The emissivity does depend on F if the den- 
sity is large enough that the gas is no longer highly ionized 



(rini/nn - 0.5«h[7'/10'^ KJ-O-^^p-j', e.g. Schaye 2001b). 

3.2. Self-shielded Gas 

If the column density of a cloud is sufficiently high, gas 
in the central regions becomes self-shielded from the meta- 
galactic background and ionizing radiation cannot penetrate 
the cloud. Unfortunately, because the simulations do not 
include radiative transfer, we cannot identify these clouds 
self-consistently. Instead we note that photoionized, self- 
gravitating clouds in (local) hydrostatic equilibrium obey 
iisd^2001bf) 

ilk) 

(1) 

where A^hi is the column density of neutral hydrogen. (Here, 
we have set the mass fraction of gas in the cloud to fib /ft,,,.) 
If we assume that the gas in the simulation is close to hydro- 
static equilibrium, we can then identify self-shielded clouds 
through the physical densities reported by the simulation. The 
clouds should become self-shielded above some column den- 
sity threshold that depends on the shape of the ionizing back- 
ground (because higher energy photons have smaller absorp- 
tion cross-sections). To fix this density, we follow 'Schay^ 
(2001a), who performed detailed radiative transfer through 
such clouds. These calculations show that self-shielding be- 
comes significant at A^hi ^ 10''^ cm^^ (see also Katzetal] 
1996a). At z = we find this column density to correspond to 
a physical density threshold «h,ss = 10^^ cm^^. We also find 
«H.s,s oc r^/-' in the range of interest, as predicted by equa- 
tion Q; the relatively modest evolution in the spectral shape 
of the Haardt & Madau 1.2001 ) ionizing background has lit- 
tle effect on the result. For simplicity, we will use this scal- 
ing to determine how nH„ss evolves with redshift. Of course, 
dense gas will not self-shield if it is collisionally ionized. To 
crudely account for this, we will assume that gas must have 
T <T,s = 10*-^ K to be self-shielded. 

Now that we have identified it, we must compute for 
shielded gas. In principle, this is straightforward. With- 
out an ionizing background, most of the gas will be in ap- 
proximate collisional ionization equilibrium (CIE), and 
is again simply a function of the local density and tempera- 
ture. In this case, we compute with Cloudy through a grid 
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over temperature; the spacing is AlogT ~ 0.02 in the range 
3.7 < logT < 4.3 and AlogT = 0.05 elsewhere (where T is 
in degrees K). The recombination rate at T < lO^*'^ K depends 
slightly on metallicity, so we also compute a grid in metallic- 
ity with spacing of one dex. However, is so small in this 
regime that the metallicity turns out to have only a negligi- 
ble effect on our results. We do not need a grid in density, 
because CIE. Note that, except at the highest tem- 

peratures, Lya following collisional excitation dominates by 
a large amount over Lya emission following recombinations. 
We show ta for gas in CIE as the dot-dashed curve in Fig- 
ure[n Note that it has a much steeper temperature dependence 
than optically thin gas has. We shall see that this has impor- 
tant consequences for our different treatments of self-shielded 
gas. 

The simplest prescription, and our first case, is then to take 
particle temperatures reported by the simulation at face value 
and assume CIE. Unfortunately, this is probably not accurate. 
As Fardal et al. (2001) point out, the introduction of an ioniz- 
ing background into any simulation without a self-consistent 
treatment of radiative transfer leads to an unphysical energy 
exchange between self-shielded gas and the ionizing back- 
ground. In our simulation, the thermal evolution of each par- 
ticle is computed assuming that it is exposed to a uniform 
ionizing background that photoheats the gas. As we shall see 
below, cool dense gas approaches an asymptotic temperature 
at which photoheating nearly balances radiative cooling. Al- 
though self-shielded particles can still cool radiatively, they 
(by definition) do not experience any heating from the ioniz- 
ing background. The simulation therefore overestimates their 
temperature. Using the model of Schave (2001a), we find 
that if photoheating from the metagalactic background were 
the only heat source, these particles would quickly cool be- 
low T ^ 10^ ' K and produce little or no Lya emission. Thus 
our second case - and the most conservative - is to set Eq, = 
for self-shielded gas. Of course, the temperature could remain 
high because of other heat sources that are not included in the 
simulation. 

On the other hand, it is possible that some fraction (or even 
most) of the gas that we label "self-shielded" is actually pho- 
toionized. Nearly all self-shielded regions in the simulation 
surround galaxies, so the local ionizing background may ex- 
ceed the mean by a large factor Our third case is thus to 
assume that all gas is optically thin and to use the simulation 
temperatures. (Note that the simulation temperatures could 
still be wrong if the local ionizing radiation differs from the 
mean background. However, the temperature dependence is 
relatively modest for optically thin gas, so such errors are not 
as important in this prescription.) Our first case, in which we 
assume CIE, is in some sense between the second and third: 
the local ionizing background is strong enough to heat the 
gas, but weak enough that CIE is a good approximation. Ob- 
viously, none of these solutions is either elegant or fully self- 
consistent, but we expect that they should bracket reality. 

We note that Fard al et aTl (1200 II) took a different approach 
and used a simulation without an ionizing background. In this 
case the particle temperatures are fixed exclusively by gravi- 
tational processes and galactic feedback. If photoheating can 
be neglected before the gas becomes self-shielded, this would 
yield more accurate temperatures for the dense, shielded gas. 
However, we will show below that photoheating is usually 
important for the Lya emission. 

Our treatment of cool, dense gas has a number of additional 
uncertainties. First, in computing its emissivity, we assume 



that the gas is optically thin to Lyman line radiation. This is 
obviously not true; in fact, many photons from the higher Ly- 
man lines will be absorbed inside the cloud and downgraded 
to either a single Lya photon or a pair of 2s continuum pho- 
tons through radiative cascades. This could increase the Lya 
emissivity by a factor of a few in the best cases. Unfortu- 
nately, the fraction of these photons that are reprocessed into 
Lya photons depends on the column depth of the cloud, its 
geometry, and its velocity structure (because photons escape 
when they scatter into the wings of the Lyman lines). Given 
the other uncertainties, we do not attempt to model this ampli- 
fication. Second, we neglect absorption by dust in the cloud; 
this is particularly important if the Lya photons scatter many 
times before escaping. Our simulations do follow the metal- 
licity of each gas particle, so we could in principle estimate 
this effect. But it too depends on the geometry and veloc- 
ity structure of the cloud (e.g., Neufeld 1991), and we will 
not attempt to do so. Because the IGM metallicity is nor- 
mally modest, we expect dust absorption to be less important 
in these environments than inside galaxies. 

3.3. Star-forming Particles 

So far we have considered gas outside of galaxies. Ac- 
tive star formation (SF) also produces substantial Lya emis- 
sion when gas in the host galaxy absorbs ionizing radi- 
ation from stars and subsequently recombines. With a 
slight abuse of terminology, we will refer to this mecha- 
nism as the SF component. In the simulations, SF occurs 
in any gas particle whose den sity exceeds n^ = 0.129 cm"-' 
( Springel & Hernauist 2003a). This threshold was fixed 
through comparison to observations of star-forming galax- 
ies (Kennicutt 1989, 1998), which show an analogous sur- 
face density threshold for star formation. Independent mod- 
els of self-gravitating galactic disks exposed to UV radiation 
find that a gravitationally unstable cold phase begins to appear 
above a fixed physical density threshold, although the precise 
value is somewhat uncertain (Schave 2004). If we decreased 
the threshold, star formation would be more widespread and 
some of the emission we attribute to self-shielded gas would 
instead come from inside of galaxies (and be dominated by 
that due to star formation). 

For each star-forming particle, we convert the SFR re- 
ported by the simulation to a Lya luminosity via La — 
10'*^ (SFR/M0yr^') erg s^', which is accurate to within a 
factor of a few according to the stellar population synthe- 
sis models of Leitherer et al. ( 1999) for SF episodes with a 
Salpeter IMF, a stellar mass range of 1-100 M^, a nd metal- 
licities between O.OSZq < Z < 2Z0 dLeitherer et al.tl999) . as- 
suming t hat ^2/3 of ioni zing photons are converted to Lya 
photons ('Osterbrock'1989). Our assumption of case-B recom- 
bination should be valid provided that most ionizing photons 
are absorbed in the dense interstellar medium of the galaxy. 
The actual Lya luminosity of a galaxy depends strongly on 
the distribution of ionizing sources, the escape fraction of ion- 
izing photons, the presence of dust, and the kinematic struc- 
ture of the gas (e.g., Kunth et al. 2003), so this should be taken 
as no more than a representative estimate. Dust is particu- 
larly important, because Lya photons must travel a long dis- 
tance before scattering out of resonance. For example, only 
~ 20-25% of bright star-forming galaxies at z = 3 have high- 
equivalent width Lya lines (SOO; Shaplev et al. 2003), with 
the rest showing weak or no emission. Because our simula- 
tions do not include this possibility, we expect our results to 
overpredict the abundance of bright Lya emission by a com- 
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parable factor. Note that we assign the entire Lya luminos- 
ity to the star-forming particle, which amounts to assuming 
that all ionizing photons are absorbed close to their source. 
In reality, some fraction will escape to distant regions of the 
galaxy or even into the IGM. Such detailed radiative transfer 
is beyond the capabilities of our simulations, so we will sim- 
ply note that Lya emission from star formation could be more 
widely distributed. In this case, although the energy source 
is inside the galaxy, the Lya photons would still allow us to 
probe the IGM. 

In principle, some fraction of the interstellar medium of 
galaxies could also produce diffuse Lya emission, especially 
if it is heated to T - lO'* ^ K. A gain, there is no way to follow 
the detailed structure of the ISM in our relatively coarse cos- 
mological simulations. We will therefore neglect such emis- 
sion processes. 

3.4. Pumping from Higher Order Lyman Lines 

The thick curves in Figure [2 neglect Lya emission from 
radiative pumping (i.e., absorption of a Lyman photon fol- 
lowed by a radiative cascade that produces a Lya photon). In 
Cloudy, we achieve this using the "no induced processes" op- 
tion. One reason for this choice is that we wish to exclude 
Lya photons emitted immediately after absorbing a Lya pho- 
ton from the background radiation field: this process does 
not contribute to the net luminosity of the gas. However, we 
should include absorption of higher Lyman line photons that 
cascade to Lya. 

Absorption of a Ly/3 photon {Is 3p) cannot result in 
emission of a Lya photon {2p Is): the Ly/3 photon puts 
the atom into the 3p state, from which it must decay to the 
li or 2i state. Unfortunately, Cloudy assumes that all levels 
with n > 3 have their orbital angular momentum states com- 
pletely mixed, w hich is only a good assumption for densities 
nil > 10** cm"^ jPengellv & Seatonlll964ft . so it does allow 
unphysical conversion of Ly/3 to Lya. This is another reason 
why we choose the "no induced processes" option. 

The selection rules do allow all higher Lyman lines to cas- 
cade to Lya. As an example we now estimate the impor- 
tance of Ly7 pumping, in which the atom begins in the 4p 
state. The selection rules allow: (1) emission of another Ly7 
photon through direct decay to the ground state, (2) emis- 
sion of two photons following decay to 2s, or (3) the chains 
4p 3s 2p Is and 4p ^ 3d ^ 2p ^ Is, which both 
produce Lya photons.^ Thus the fraction of absorptions re- 
sulting in a Lya photon is 

f _ A4p^2d +A4p3s 

J-f.a — . . . . ) (■^) 

^4p,3d +^4p,3s +A4p2s +^4^,15 

where A v,,. is the Einstein A-coefficient between states x and y. 
The result is „ — 0.0444, where we have averaged over the 
4p states. This is the fraction for a single absorption; however, 
if the medium is optically thick to Ly7 photons, a re-emitted 
Ly7 photon will be repeatedly absorbed until it is transformed 
into a Lya photon or a pair of 2s continuum photons. To ac- 
count for this possibility, we could set A^p is ^ in equation 
which yields 0.278 L 

The thin lines in Figure show the extra Lya emissivity 
from Ly7 absorption for several different densities at z = 3, 
assuming an optically thin medium. We see that in this case 
pumping is only a small correction for all densities; on the 

^ Note however that not all 4p states allow decay to 3d because the total 
angular momentum can only change by ti. 



Other hand, if the medium were optically thick to Lyman 
line photons, the correction could be comparable to the emis- 
sivity from recombinations in the low-density IGM.^ How- 
ever, it is always much less than the emission from colli- 
sional excitation in dense gas with T > 10*-^ K, which, as 
we will see, dominates the total. For simplicity, we will thus 
neglect this process in the following and note only that it 
could increase the luminosity of the low-density and/or low- 
temperature IGM by a small factor. Including even higher 
transitions than Ly7 does not change this general conclusion. 
However, it should be noted that pumping could be much 
more important in places where the intensity in the Lyman 
lines significantly exceeds that of the mean background radi- 
ation. 

Finally, we note that the pumping contribution increases 
with redshift, because in photoionization equilibrium the 
mean neutral fraction increases with the proper density. 

4. IGM CHARACTERISTICS 

In this section we examine some general properties of Lya 
emission i n the simulations, including the phase diagram of 
the IGM (^4.H . the underlying energy source for the Lya ra- 
diation ( ^4-.2\ . and the effects of galactic winds (' ^4.3> . 

4.1. Phase Diagram of the IGM 

Figure |2] shows nu-T phase diagrams for two simulation 
outputs: G6 at z = (top row) and Q5 at z = 3 (bottom 
row). The left column shows the particle number density 
(unweighted), while the right column weights by the Lya 
luminosity of the particles, assuming CIE for self-shielded 
gas. F03 describe most of the main features in detail. The 
majority of particles lie along a line connecting cool, un- 
derdense gas and moderately overdense, warm gas. This is 
the di ffuse, pho toionized IGM responsib le for the Lya for- 
est (K atz etalJl 9 %: Hui & Gnedin 1997l lSchave etal Jl99qt: 
iMcDonald et alJ 2001). Another set of particles occupy a 
broad range of temperature at moderate to large overdensi- 
ties. This gas is the shocked IGM (now known as the "warm- 
hot IGM" at z ~ 0). At z = 3, much of this gas has been 
shocked by galactic winds, but by z = the nonlinear mass 
scale has risen significantly and large-scale structure shocks 
provi de most of the heating (Cen & Ostriker 1999; Dave et aD 
2001: Croft etal. 2002; K eshet et al. 2003; Furlanetto & Loebl 
.2004J . A third set of particles lies along a nearly vertical line 
at r ~ 10^ K, representing gas that has cooled and collapsed 
into bound objects but is not currently forming stars. We will 
refer to this feature as the "cooling locus." In simulations 
without winds, the vast majority of these particles have not 
yet formed stars but will do so in the relatively near future. 
With winds, the interpretation is more subtle (see ^4.31 be- 
low). The rapid decrease in the cooling rate at T < 10"*^ K 
forces this gas to approach a constant asymptotic temperature 
near that value. The high-temperature envelope of this locus 
appears because the cooling time falls rapidly as the tempera- 
ture increases. Although very few particles inhabit this region 
of phase space, their emissivity is so large that they do appear 
in panel (b). 

It is clear that most of the Lya emission comes from gas 
on the cooling locus, especially at z = 3, because this curve 
lies near the peak of the Lya emissivity in CIE (note that the 

* To accurately compute the pumping contribution for a gas cloud that is 
optically thick in the Lyman lines we would need to do a full radiative transfer 
calculation because the Lyman line intensity would vary within the cloud. 
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Fig. 2. — Phase diagrams of the IGM. (a): Particle number density in the G6 simulation at z = 0. The dotted line shows the mean density at this epoch, (b): 
Pailicle number density weighted by Lya emission. The solid lines mark the self-shielded region of phase space; we assume CIE for gas in this region, (c), (d): 
Same as the the top row, except for the Q5 simulation at z = 3. All contours are equally spaced (logarithmically) in density, with an interval of 0.5 dex. 



majority of this region is in the self-shielded regime). Be- 
cause is so sensitive to temperature at T ^ 10'*'^ K, the 
emission characteristics will depend on the precise tempera- 
ture along this curve. This is particularly evident in Figure 
IIJ?: the densest gas has T < 10^ ' K, below the peak of Lya 
emission, and so contributes only a small fraction of the to- 
tal luminosity. In contrast, at z = 3 the temperature is higher 
and all of the gas contributes significantly. As emphasized in 
83.21 much of this cool gas will actually be self-shielded from 
the mean metagalactic ionizing background (which is respon- 
sible for most of the heating in the simulation) and hence may 
cool to a lower temperature than the simulation allows; in this 
case the emissivity could change dramatically. The sensitive 
dependence of the emission on the temperature is a direct con- 
sequence of assuming CIE. If all of the gas remains optically 
thin, is relatively flat for T < lO'*^ K (see Fig.^, so the 
densest gas at z = would still emit strongly. 

We do not show star-forming particles in Figure|2] Particles 
in the simulation form stars when hh > 0.129 cm^^, so they 
would all lie above the upper edge of these plots. 

4.2. The Energy Source ofLya Emission 

'Ha iman et all ( 2000') and ' Fardal et al.1 (|200D have sug- 
gested that gravitational shock-heating can yield large Lya 
lumino sities for halos in the midst of collapse. Haiman et al. 

argued that the gas would be shocked to the virial tem- 
perature of the halo and cool rapidly to T ~ 10"* K through 
a combination of free-free and line emission (but primarily 
higher excitation lines than hydrogen Lya). The accreted 
baryons would then be left out of hydrostatic equihbrium and 



collapse to the center of the halo. During this stage, the gas 
remains at T ^ 10"* K and radiates its gravitational energy in 
Lya photons. iFardal et al.l (120011) found that most of the gas 
was not shocked to high temperatures in their simulation, so 
that a potentially even larger fraction of t he gravitational en- 
ergy could be radiated in Lya photons. 'Birnboim & Dekell 
(2003) reach similar conclusions from analytic arguments 
and one-dimensional simulations, and K eres et al.l (2004) con- 
firmed the conclusion in a detailed analysis of gas accretion 
in three-dimensional cosmological simulations. In either case, 
we expect a large fraction of the emission to come from gas 
that has recently been accreted, and naively that gravitational 
processes provide the energy reservoir for the Lya emission. 

We can test this hypothesis with our simulations. The 
primary additional ingredient we have added is the ionizing 
background. For each particle in the simulation, we com- 
pute the photoheating rate as well as the total cooling rate, 
including radiative cooling and adiabatic expansion. The ra- 
tio of these two quantities yields the fraction of Lya energy 
that can be directly attributed to the ionizing background; the 
remaining energy can come from gravitational processes (or 
any other, such as galactic feedback). We show the results in 
Figure |3] The upper solid curves show the ratio of (radiative 
heating/total cooling) along fixed density slices. The dashed 
curves show the particle distribution at that density (with arbi- 
trary normalization). We have selected densities on the cool- 
ing locus; note that «h,ss = 10~^'^^, lO^-' cm^-' at z = 3, 0. 
The (heat/cool) ratio behaves erratically at low temperature 
because there are few particles in this regime, many of which 
have unusual histories (such as recent ejection by winds). The 
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Fig. 3. — Ratio of photoheating to total cooling along fixed density slices (solid curves). The dashed curves show the particle distribution function at this 
density (with arbitrary normalization). The vertical dotted line marks the point where the distribution function peaks; the horizontal line shows where heating 
and cooling just balance. 



distribution peaks where photoheating nearly cancels cool- 
ing. The particles are still cooling (albeit much more slowly 
because of the ionizing background), especially in the low- 
density slice, and heading toward the star formation threshold. 
Clearly the particle distribution is determined by the balance 
of photoheating and cooling. Note that this is true even for 
self-shielded gas; as we have argued in ^3.21 such behavior is 
an important limitation of the simulations, as it represents an 
unphysical heat supply for such gas. 

The solid line in Figure 0] shows the fraction of Lya emis- 
sion due to photoheating when we integrate over the entire 
particle distribution (assuming Ea — for self-shielded gas). 
At ^ ^ 1' photoheating provides somewhat more power than 
gravity or winds (accounting for ^ 2/3 of the total emission). 
On the other hand, at z ~ 0, photoheating is negligible be- 
cause the ionizing background has fallen substantially. These 
results imply that the photoionizing background is important 
for optically thin gas; however, with our fiducial Wss, such gas 
accounts for only ~ 1-10% of the total emission. Whether 
photoionizations are important for the total emission depends 
on our assumptions about cool dense gas. 

Unfortunately, it is difficult to self-consistently vary the ion- 
izing background, because it is built into the simulations. In 
particular, the equilibrium temperatures of cool dense gas par- 
ticles depend on the background. A weaker background de- 
creases the temperature and hence the emissivity; a stronger 
background will usually increase the emissivity unless it 
pushes the temperature over the peak. To include these ef- 
fects rigorously would require radiative transfer to be part of 
the simulation, which is not yet feasible. We will therefore 



take an approximate approach to varying the background. 

First, we can vary our criterion for "self-shielded" gas. This 
tells us how important the ionizing background is as a func- 
tion of density for gas on the cooling locus. It also mimics 
a situation in which the local ionizing radiation field is larger 
than the mean. The error bars in Figure |3 show what hap- 
pens if we increase or decrease «ss by 0.67 dex. Because the 
cool dense gas has its temperature fixed almost entirely by the 
ionizing background, increasing the threshold density also in- 
creases the fraction of the emission from this mechanism. The 
dot-dashed line shows the fraction if n^^ — 0.129 cm^^, the 
SF threshold; in this case, > 90% of the Lya emission comes 
from the ionizing background. We conclude that photoioniza- 
tion will strongly dominate gravitational heating as the source 
of Lya emission, regardless of redshift, if the local ionizing 
field is strong. 

Second, we can change the amplitude of the ionizing back- 
ground when calculating the emissivity and heat/cool ratio 
(but not the gas temperature, because that is built into the sim- 
ulation). The dashed curves show what happens if we increase 
or decrease T by one order of magnitude. This choice also 
modifies n^s by 0.67 dex (see equation Q) as well as modi- 
fying the distribution of emissivity by decreasing the impor- 
tance of recombinations following ionizations. In particular, 
while strengthening the ionizing background appears to have 
little effect other than on Wss (except at z 5), weakening it re- 
duces the importance of the UV background by a much larger 
amount. This is because with such a background most of the 
emission arises from hot, dense gas at T > lO"*'^ K for which 
collisional excitation dominates completely (see Fig.Q- 



8 




12 3 4 



z 

Fig. 4. — Fraction of Lya radiation ultimately due to the ionizing back- 
ground. The solid curve shows the fraction for gas below our fiducial riss; 
error bars show the fraction if we increase or decrease the threshold by 0.67 
dex. The dot-dashed curve assumes that all gas is optically thin. The upper 
and lower dashed curves increase and decrease the ionizing background by 
one order of magnitude from our fiducial choice. 



With these resuhs in mind, we now consider four cases for 
emission from cool, dense gas. First, it could remain optically 
thin because of photoionization by embedded sources. Then 
the simulation is adequate: the temperatures are incorrect be- 
cause the amplitude and shape of the ionizing field change, 
but the temperature is relatively unimportant for optically thin 
gas. In this case, the dot-dashed curve in Figure |4] shows that 
photoionization completely dominates. Second, self-shielded 
gas could be invisible, either because it cools rapidly or be- 
cause of dust. In this case the simulation is exact (because 
only optically thin gas is observable), and the photoionizing 
background dominates at z > 1 although not to the exclusion 
of gravitational processes. Third, self-shielded gas remains 
visible and is heated by some unspecified source (this cor- 
responds to our CIE curves). The simulation temperatures 
are certainly unreliable in this case. While optically thin gas 
would still be powered by the ionizing background, the self- 
shielded gas (and hence most of the emission) would by defi- 
nition be powered by some other mechanism. Fourth, the gas 
could cool once it condenses beyond Wss and remain visible 
only during the initial cooling. Here, we have found that the 
ionizing field plays an important (though not necessarily dom- 
inant) role in the initial temperature (except at z = 0). Deter- 
mining which of these possibilities corresponds to reality re- 
quires more advanced simulations that include self-consistent 
radiative transfer. 

This does not imply that shocks are irrelevant to structure 
formation: as 'Haim an et al.l (l200Cft point out, a large fraction 
of the shock energy could be radiated through other transi- 
tions. We briefly consider the importance of shocks here. We 
have randomly selected luminous particles from t he z — 2.75 
snapshot of the 03 simulation. (As we will see in ^4.31 winds 
complicate the interpretation of particles on the cooling lo- 
cus, so we have here chosen our highest resolution simula- 
tion without them.) For this simulation, we have snapshots 
spaced at Az = 0.05 for z = 2.75-3.25 and at Az = 0.5 for 
3.5 < z < 7. For each particle, we find the largest temperature 



Tmax and entropy K^nax (where K ~T / rp-/^) it has experienced 
in this set of snapshots. (We exclude any snapshot at which 
the particle formed stars, but these are extremely rare any- 
way.) Note that this is actually only a lower limit to the true 
maximum, because the outputs are spaced fairly coarsely and 
the coolingtime is short for large temperatures and densities. 

Figure js] compares the maximum and current T and 
K for particles with —27.25 < loge^ < —26.75 (trian- 
gles) and —25.25 < logec < —24.75 (crosses, in units 
of erg cm^^ s^'). The former set has a relatively low emis- 
sivity and is not self-shielded; the latter is well within the 
self-shielded regime. We find that most of the lower emis- 
sivity particles have Tmax > r(z = 2.75) (some much greater 
than this level). We expect virial temperatures of ~ 
105(M/10"' M0)2/3 Katz = 3 (Ba rkana & L oeb 2001): most 
of the resolved star-forming halos in this simulation have 
masses above 10^° Mq. Interestingly, most of the particles 
remain at or belo w 2.5 x 10^ K, the threshold defined by 
iKeres et alJ ( l200l to divide accreting gas into a "cold" and 
"hot" mode. They argue that virial shocks are responsible for 
the latter phase, while cold accretion occurs along filaments 
and is subject only to weak shocks. Our conclusions are con- 
sistent with theirs; however, many of our halos have virial 
temperatures compar able to the threshold, so the division is 
not as clean. Like Keres et alJ (2004), we do see some evi- 
dence for "bimodality" in the temperature distribution, with 
very few particles having Tmax ^ 10^^ K. This is likely be- 
cause the cooling function has a peak at these temperatures. 
Typically, particles in this emissivity range also have /Tmax 
several times larger than K{z — 2.15), indicating a significant 
amount of radiative cooling regardless of the accretion phase. 

The higher luminosity particles show smaller maximum 
temperatures: most have T^ax ^ lO'^'^-lO^ K, but many have 
remained even cooler. They have also lost more entropy, from 
initial levels comparable to the other emissivity bin. This is 
indicative of substantial radiative losses. One possible expla- 
nation for the difference is that many of these particles were 
heated at z > 3.25, so that we are more likely to have missed 
the shock in our snapshots. Alternatively, these may have 
been accreted through the cold mode identified by Kere s et all 
(2004). Those authors argue that the cold mode is ultimately 
responsible for most of the star formation at z ^ 3, and our 
results appear consistent with this conclusion. We empha- 
size again that, because of the short cooling times and limited 
number of outputs, we actually only have lower limits on Tmax 
and/iTmax- 

In summary, we cannot neglect the ionizing background 
(and any additional contribution from local sources) when 
computing the Lya emission from the IGM. Unfortunately, 
its precise role depends on self-shielding and local radia- 
tion sources, so a self-consistent treatment will require ra- 
diative transfer in full cosmo logica l simulations , including 
its dynamical effects. Haima n et alJ ( 12000) and .Fardal et alJ 
(|2p01 ) neglected the ionizing background and argued that the 
gas would contract after being shocked and radiate its grav- 
itational energy. Instead, we find that the gas settles into 
a quasistatic phase where photoheating nearly balances ra- 
diative cooling. This heat exchange becomes unphysical if 
(at some point during the contraction) the gas becomes self- 
shielded (see ^3.2i . The density at which this transition oc- 
curs will determine how much energy the metagalactic ion- 
izing background supplies compared to gravitational collapse 
(and winds). 
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Fig. 5. — Left: Compai'es the cun'ent temperature of simulation pailicles to the maximum temperature achieved in the previous simulation outputs. Since we 
only have a finite number of outputs, the plotted maximum temperatures are lower limits to the true maxima. We have randomly selected 100 particles from the 
03 simulation at z = 2.75 for each emissivity bin; the triangles and plus signs have —27.25 < \ogta < —26.75, and —25.25 < \ogta < —24.75, respectively 
(in units of erg cm~' s~ ' ). The latter bin is composed of self-shielded pailicles (for which we assume CIE). Right: Same, but for the entropy. 



4.3. Winds 

Most of the simulations th at we study include a galac- 
tic wind mechanism (see Soringel & H ernquistI l2003al: 
iFurlanett'o'et al. 200 4a for detailed descriptions of the wind 
implementation). We might expect this prescription to have 
important consequences for the Lya emission. We have there- 
fore also examined some lower resolution runs with weak or 
no winds (P3 and 03, respectively). In this section we review 
the effects of the wind prescription on the simulation particles. 

The most important consequence of winds is to sharply 
reduce the star formation rate: the SFR density is 
(1.18,O.855,O.244)/i3M0yr"i Mpc"^ at z 3 in the (03, 
P3, Q3) simulations, bringing the Q3 sim ulation into rea- 
sonable agr eement with observations ( Springel & HernquistI 
l2003bt iHer nauist & Springel 2003). The winds accomplish 
this in two ways. First, in the absence of ongoing accretion, 
they simply deplete the gas in the star-forming phase by eject- 
ing it into the IGM. Second, winds interact with gas currently 
accreting onto halos. In small galaxies, a wind can severely 
reduce the accretion rates by entraining gas and transporting it 
out of the halo or by shedding its energy to gas in the outskirts 
of the halo. In large galaxies, on the other hand, winds remain 
bound to the halo and only weakly affect the accretion rate. 
The halo gas rids itself of any wind energy through cooling. 
In this case even the wind material itself can re-enter the star 
forming phase after a short delay. 

Thus, winds have a substantial effect on gas lying along 
the cooling locus. In their absence, nearly all of this gas is 
pristine (i.e., has not formed stars in the past). If we trace 
the evolution of individual particles, they slowly cool as their 
densities increase until they are able to form stars. If we turn 
on weak winds (i.e., those that are unable to escape the grav- 
itational potential wells of typical galaxies), we would expect 
more gas to lie on the cooling locus because winds eject star- 
forming particles into the halo and provide extra heat to the 
accreting gas. Moreover, a significant fraction of the dense 



gas would have formed stars in the past and thus be metal- 
enriched. As the winds strengthen, the amount of cool gas 
should decrease again once they can completely eject parti- 
cles from galaxy halos. However, some fraction of the wind 
particles will still be trapped along "accretion channels" and 
fall back on the galaxies. 

We find precisely this trend in the simulations: the 03 sim- 
ulation actually has the least amount of gas along the cooling 
locus, even though it has the highest SFR. The Q3 simulation 
has slightly more gas in this phase, while the P3 simulation 
has the most (apparently these weak winds get trapped near 
the galaxies). In both the Q3 and P3 simulations, a signifi- 
cant fraction of high density, low temperature particles have 
formed stars in the relatively recent past. Thus winds have 
complex, and not always intuitive, implications for Lya emis- 
sion. However, the effect on the fraction of gas lying along 
the cooling locus is modest, varying by only a small amount. 
The upper left part of Figure|2l: appears essentially unchanged 
between the different wind prescriptions. 

A more obvious effect is on the shocked, diffuse IGM at 
high redshifts: without winds, this phase would be consid- 
erably weaker in Figure Of course, current cosmologi- 
cal simulations cannot fully resolve the internal structure of 
galactic winds. In the nearby universe, winds are observed to 
have complex multiphase media, including a hot diffuse in- 
terior, dense warm clouds, and (in some cases) a cool dense 
shell of entrained matter. These components are largely un- 
resolved in our simulations and could strongly increase the 
Lya luminosity of the wind region and render its geometry 
more accessible to observations. The true effect of winds may 
thus be more complex than simply adding gas to the hot, low- 
density phase. 

Winds are clearly complex phenomena, and their effects on 
Lya emission are difficult to predict a priori. It is obvious 
from Figure I2J: that the additional hot gas will make a negli- 
gible difference to the emission. We will exa mine the effects 
on bright emission from the cooling locus in ^5.3l below. 
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Fig. 6. — Maps of Lyct emission from the IGM {a and c), assuming self-sliielded gas is in CIE, and from star formation {b and d). All are from the Q5 simulation 
at 2 = 3. Sources identified by SExtractor (with the thi'eshold at $ = 100 photons s~' cm~^ sr~ ' or 9.6 X 10^^' erg s~ ' cm~^ arcsec"-) are outlined in green, 
with the size proportional to the source's area on the sky. The slice shown in panels (a) and (b) is unusually rich; the other is more typical. The field is 
2.75/;^' Mpc wide and 0.67/;^' Mpc thick (comoving; this con'esponds to 1.2 A). The angular resolution is 1". 



5. LYaMAPS 

We will now describe the characteristics of narrowband 
Lya maps in our simulations. In this section we will focus on 
qualitative features and global statistics of the maps. We will 
discuss the statistical properties of well-localized Lya emit- 
ters in |6] Figure|6]shows two example slices of the universe 
at z = 3 (taken from the Q5 simulation). The slice is 1.2 A 
thick (corresponding to Az = 10^^ or 0.67/i^' Mpc comov- 
ing), and the field of view is 2.1' with angular resolution 1" 
(these correspond to 2.75/i^' Mpc and 21.6/;^' kpc, both co- 
moving). Unless otherwise specified, all of our z = 3 maps 
assume these parameters. Panels (a) and (c) show the emis- 
sion from IGM gas (assuming that the self-shielded gas is in 
CIE), while (b) and (d) show that from star-forming particles 
in the corresponding slices. The green circles outline distinct 
sources (see ^6. 21 for details). The slice at right is typical for 
this redshift; the left slice is among the richest few percent in 
the simulation. 

These maps have a number of generic features. First, bright 
Lya emission from the IGM and from SF are highly corre- 
lated. It is rare for one to appear without the other; as we 
have described above, this is because most of the luminous 
IGM particles are either cooling slowly to form stars or have 
recently felt the effects of galactic winds. The total luminos- 
ity of each region is usually dominated by star-forming par- 
ticles. However, the IGM emission tends to be distributed 
on somewhat larger scales than the SF. Most of the emission 
from these coronae (especially in their cores) comes from 
gas above the self-shielded threshold. As we will show be- 
low, these cores depend sensitively on our assumptions about 
self-shielded gas. Second, Lya emission is aligned along the 
cosmic web of filaments, with low surface brightness emis- 
sion (from optically thin gas) connecting the coronae of bright 
emission surrounding galaxies. These properties remain true 
at all redshifts, though the importance of the coronae varies 
strongly, as we shall show next. 

5.1. Redshift Evolution 

Figure shows how the statistical properties of the maps 
evolve with cosmic time. For each redshift, we show the prob- 
ability distribution function dn /d In $ of pixel surface bright- 
ness $ computed from 120 slices of the highest resolution 



simulation available at that redshift. In order to better inter- 
pret the physical processes as they evolve, we have held the 
comoving spatial resolution and slice thickness constant with 
redshift; note then that the angular resolution and spectral 
width of the observations differ. (They equal the choices of 
Figure|6lat z = 3.) We have assumed CIE for the self-shielded 
gas, although the choice makes no visible difference to this 
Figure. At all redshifts, dn/dln$ peaks at a relatively small 
surface brightness characteristic of the median cosmic den- 
sity. In the highly-ionized low-density IGM, recombinations 
dominate over collisional excitation (see Fig.flJ, so the mean 
surface brightness can be estimated from the recombination 
rate (Hogan & Weymann 1987; Gould & Weinberg 1996). It 
is: 
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where rj w 0.42 is the fraction of recombinations that produce 
a Lya photon, Wrec is the recombination rate, h{z) — H{z) / Hq, 
and 5 represents the effective overdensity (i.e., the density rel- 
ative to the cosmic mean) of the observed pixel. Here we have 
used a = 4.2 x 10^ '"^ cm-' s^' for the (case-A) recombina- 
tion rate, appropriate for gas with T — 10"* K. (This estimate 
does not apply to the bright coronae, which are not necessar- 
ily highly ionized and which have broken off from the cosmic 
expansion.) We see that $ increases with redshift because the 
mean density and recombination rate increase rapidly. This 
equation yields a reasonable estimate of the peak of dn/d<I>, 
although note that the appropriate 6 changes as structure for- 
mation progresses. 

Unfortunately, the probability distributions of the pixel flux 
peak many orders of magnitude below realistic detection 
thresholds, regardless of redshift. We therefore focus on the 
bright end of the distribution in Figure |8] where we show the 
fraction of pixels above a given threshold surface brightness, 
F{> $). The first three panels show emission from the IGM 
if self-shielded gas is in CIE, if it has = 0, and if it is op- 
tically thin to the ionizing background. The fourth shows the 
total emission from the map (including SF, and assuming CIE 
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for the self-shielded gas). 

Looking first at panel (a), we see that the bright pixels 
evolve in a way that differs qualitatively from the low-density 
IGM. The brightest emission peaks at z 2-3 and decreases 
rapidly toward high and low redshifts. The decline occurs 
when the cooling locus shifts to low temperatures because the 
ionizing background is weak and/or soft. As a consequence, 
declines rapidly if we assume CIE, so the maximum pixel 
brightness falls. The effect is particularly severe at z ^ 5, 
where the temperature of this component is T « 10'* '^ K and 
self-shielded gas is essentially invisible. 

Comparison to panel (b), where we set — for self- 
shielded gas, shows that nearly all of the emission with $ > 
100 photons s^' cm^^ sr^' comes from self-shielded gas. If 
this component cannot emit (either because of strong dust ab- 
sorption or because the simulation overestimates its tempera- 
ture), the Lya coronae become orders of magnitude weaker. 
Clearly the choice of n^s is crucial in the no emission case. If 
our prescription is pessimistic (because of embedded sources, 
for example), the coronae rapidly strengthen. F03 show an 
explicit example of how the z = statistics change with the 
self-shielding cut. Thus observations of diffuse Lya emission 
around galaxies can strongly constrain the behavior of cool 
gas in these environments. 

Figure Isj:, on the other hand, shows the emission if we as- 
sume that all the gas is o ptically thin to ionizing radiation. 
Although ISchavd fcoOlalh) has shown that self-gravitating 
gas clouds become self-shielded from the metagalactic back- 
ground at these densities, most of these particles surround 
star-forming galaxies. As a result, the local ionizing field 
could be much stronger than the mean, in which case the 
gas may remain optically thin. In this scenario, the maxi- 
mum brightness decreases, because optically thin gas is al- 
ways highly ionized and the CIE collisional excitation peak 
disappears. However, the fraction of moderately bright pixels 
remains nearly the same, because the same particles emit (and 
remain reasonably bright). The exception is at z = 5, where 
the optically thin case has much stronger emission. This is be- 
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diation, (d): Total Lya emission from both the IGM (including self-shielded 
gas as in panel a) and SF. 



cause the temperature dependence of is much weaker in the 
optically thin case, so our predictions are much less sensitive 
to the location of the cooling locus. 

The final panel shows that emission due to SF dominates 
the brightest end of dn/d$ at all redshifts. This is not sur- 
prising, as SF produces Lya photons efficiently. We find 
that bright emission from SF rises to z ~ 3 and then de- 
clines slowly at higher redshifts. Although the SFR slowly 
increases to z 6 in these simulations ( Scringel & Hernauisa 
l2003btlHernauist & SpringeL2003.) . the evolution flattens be- 
yond z ~ 3 because the increasing cosmic distance causes the 
maximum surface brightness to decline. On the other hand, 
note that the fraction of moderately bright pixels does not 
change significantly between panels (a) and (d) (except at 
z — 5). This is because the IGM Lya emission is more ex- 
tended (and fainter) than the star-forming regions, so IGM 
emission dominates the sky coverage at all but the highest 
surface brightnesses. Of course, we have assigned all of the 
Lya photons from SF to the host particles. Some fraction 
will actually escape, spreading this component over a wider 
area. This mechanism will help to boost the IGM emission, 
even if the ultimate energy source is (photoionization from) 
SF. Moreover, we have neglected dust absorption, which will 
likely remove a large fraction (> 75% according to SOO) of 
the SF Lya photons. Thus the bright tail may not be as pro- 
nounced as indicated in Figure |8ji and IGM emission (which 
probably suffers less dust extinction) may be even more im- 
portant. 

Figure |9] shows the redshift evolution in a different way. 
Here we plot the volume-averaged emissivity in the Lya line. 
The solid line shows the total (including both SF and the 
IGM, assuming CIE for self-shielded gas). The short- and 
long-dashed curves show £„ from the IGM if we assume CIE 
for self-shielded gas and that all gas is optically thin, respec- 
tively. The dotted curve again shows the total emissivity, but 
this time in comoving units. We have used the highest res- 
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Fig. 9. — The evolution of the volume-averaged emissivity Ea in physi- 
cal units. The solid line shows the total emissivity (assuming CIE for self- 
shielded gas). The short- and long-dashed curves show from the IGM 
if we assume CIE for self-shielded gas and that all gas is optically thin, re- 
spectively. The dotted line shows the total emissivity per comoving volume. 
Triangles show the emissivity estimated from lower resolution simulations, 
while hexagons show it with different wind prescriptions. Filled and open 
symbols refer to total and IGM emission, respectively. 



DO 



1 1 1 1 1 1 1 1 1 1 M 1 1 1 1 1 1 l_ 

\ (a) ; 

^^^^~:^r>^^^ IGM - 

~ Res \\ ~ 
-, , , 1 , , , 1 , , , 1 A, 


J 1 1 1 1 M M 1 1 1 1 1 1 1 1 1 1 

; (b) : 

SF - 

Res 

■ , , , 1 , , , 1 , , , 1 , , 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i_ 

1 (c) ; 

?='^*=»^5,^^^^ IGM - 

~ Winds V " 
" A i r ■ 


: (d) ; 

SF - 

- \ \ ^ - 
\ N 

Winds \ - 

■ r 



1 2 3 4 1 2 3 4 5 
log $ (photons cm"^ s"' sr"') 
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short-dashed curves have strong, weak, and no winds, respectively, (d): Same 
as (c), except for emission from SF. 



olution simulation available at each redshift to construct the 
curve; the triangles show estimates from lower resolution sim- 
ulations with the same parameters. The predictions have con- 
verged nicely. 

One of the crucial lessons of Figure |9lis the difference be- 
tween the two dashed curves. Evidently self-shielded regions 
dominate the IGM emissivity. The optically thin emissivity 
closely traces the total (and hence the SFR). Again, this is be- 
cause depends only weakly on temperature for optically 
thin gas, so the coronal emission depends (to zeroth order) 
only on the rate at which galaxies accrete gas (or eject gas 
through winds), which is itself roughly proportional to the 
SFR. In contrast, in CIE ea depends sensitively on tempera- 
ture. At z 2-3, quasars make the ionizing background rather 
hard and intense, so the cooling locus shifts to high tempera- 
ture and the CIE emissivities are large. But at z = and z = 5, 
the temperatures are small and CIE emissivities fall below the 
optically thin values. We will focus on z = 3 in much of the 
following, and we urge the reader to keep in mind that the CIE 
assumption is the most optimistic at this epoch. 

5.2. Simulation Resolution 

We now consider the degree to which the finite simula- 
tion resolution affects our results. Figures [Toh andFTOb show 
F(> $) for five different simulations at z = 3. The two panels 
show results for IGM emission (assuming self-shielded gas is 
in CIE; other cases yield similar results) and that due to SF, 
respectively. All the parameters of the slices are the same as 
in Figure |S] We also use the same number of slices for each 
simulation. Thus we only sample a small fraction of the vol- 
ume of the D5 and G6 boxes, and we avoid contamination 
from rare objects that do not appear in the smaller boxes. The 
figure shows that our two highest resolution runs, Q4 and Q5, 
appear to have converged. However, the lower resolution runs 
(especially G6) show an excess of extremely bright pixels to- 



gether with a deficit of moderately bright pixels, for both the 
IGM and star formation. T his may seem surprising, because 
'Springel & Hernquist'('20b3b) showed that the global SFR has 
converged between all of these simulations. We find deeper 
convergence problems because we require the spatial distri- 
bution of material within halos. In the lower resolution boxes, 
small halos are not as well-resolved and have SF (and cool 
gas) confined to only a few particles that fill more compact 
central regions. Thus the lower resolution simulations can 
have "cuspier" surface brightness distributions that underesti- 
mate the physical extents of the bright Lya regions. The lower 
resolution simulations also miss small halos, of course. On 
the other hand, it is possible that the extra large-scale power 
in the D5 and G6 simulations genuinely increases the num- 
ber of bright pixels, by (for example) concentrating the star 
formation in smaller, highly-biased regions. 

5.3. Winds 

As described in ^4.31 winds have a substantial impact on the 
SF properties of the simulation, and we would naively expect 
them to have comparable effects on the Lya emission. Fig- 
ure shows that winds do indeed have substantial impli- 
cations for the Lya emission from star-forming particles: the 
simulations with weak and no winds have nearly an order of 
magnitude more pixels with $ > lO'* photons s^' cm^^ sr~' 
than the simulation with strong winds. However, Figure [Tob 
shows that the effects on the IGM emission are quite modest. 
(Here we have assumed CIE for the self-shielded gas, but the 
differences are small for the other treatments of self-shielded 
gas as well.) The simulations with strong and no winds have 
nearly the same F{> <&), while weak winds strengthen the 
emission by only a small amou nt. T his follows the trends ex- 
pected from the discussion in ^4.31 in which we argued that 
weak winds have the strongest effect on the cooling locus be- 
cause the winds cannot escape the potential wells of their host 
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galaxies. It is not, however, obvious why the net effect should 
be so small. It is unclear whether there is a deeper reason for 
this behavior We do note that the nature of the bright parti- 
cles changes between the three simulations. Without winds, 
the vast majority of the emission comes from particles that 
have not yet formed stars and contain no metals. With strong 
winds, fewer than 50% of bright particles are pristine, because 
a number of formerly star-forming particles have mixed with 
the gas in the cooling locus. 

We also note again that the simulations do not resolve the 
internal structure of the winds. If the wind medium frag- 
ments into dense clumps inside of a more rarefied interclump 
medium, the total amount of emission could change substan- 
tially. The clumps would initially have larger emissivity be- 
cause of their increased recombination rate; as a result, they 
would cool more rapidly and could eventually pass below the 
Lya excitation threshold. As with self-shielded gas, the re- 
sulting emissivity depends sensitively on the equilibrium tem- 
perature, which would in turn depend on the ionizing back- 
ground, continued heat input from the wind, and conduction 
from the hot intercloud medium. More detailed simulations 
of individual winds are necessary to resolve these processes. 

Figure |5] also compares the volume-averaged across 
wind simulations. The hexagons show the mean emissivity for 
the P3 (weak wind) and 03 (no wind) simulations. While the 
IGM emission is approximately constant in the three cases, 
the mean SFR rises by a factor of 4-5 from the strong wind 
case, increasing the total ea by a large amount. 

5.4. Line Widths 

We now examine the velocity widths of these Lya lines. To 
estimate the distribution, we construct simulation maps as be- 
fore. For each particle, we also record its velocity (including 
both the Hubble flow and the peculiar velocity) and velocity 
dispersion (assuming pure thermal broadening). From these 
we compute the flux-weighted velocity dispersion cr (i.e., the 
standard deviation of the velocity distribution) of each pixel 
in the maps, using the same pixel sizes as in Figure |6] Note 
that a includes thermal broadening, peculiar velocities, and 
the Hubble flow between particles (but not the Hubble flow 
gradient within individual particles, because that is negligible 
for bright, compact particles). 

Figure ^2 shows the distribution of a per pixel (including 
IGM emission only, and assuming CIE for the self-shielded 
component) as a function of surface brightness $. Note that 
the slice has a velocity thickness of « 300 km s^'; uniform 
Hubble flow would yield cr « 80 km s^'. This is indeed the 
peak of the distribution in the low-density IGM. 

We see that bright pixels are narrow. Essentially all pixels 
brighter than lOphotonss^' cm^- sr^' havecr< 100 km s^'. 
At modest brightness (~ 100 photons s^ 
distribution peaks at cr « 10 to 15 km s 
accounted for by thermal broadening at the temperature of 
the cooling locus. As the brightness increases, the distribu- 
tion broadens to somewhat larger velocity dispersions, indi- 
cating that peculiar motions become relevant. This is prob- 
ably because the most luminous objects are associated with 
more massive halos with large internal velocity dispersions 
(see ^6.1> . Note, however, that the simulation has only lim- 
ited resolution for the internal structure of halos, so there 
may be extra broadening due to rotation, disc formation, etc. 
Moreover, the coronae should be optically thick to Lya pho- 
tons. In this case the photons can only escape by scatter- 
ing into the wings of the line. Thus the velocities reported 
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here will underestimate the true linewidth. For some geome- 
tries, radiative transfer can even change the line's shape (e.g., 
IZheng & Miralda-Escudell2002l and references therein). The 
best we can say is that the source regions are intrinsically nar- 
row before the effects of radiative transfer have been included. 



5.5. Comparison with the Diffuse Backgrounds 

The observability of Lya emission, especially at the low 
surface brightness levels relevant for most of the IGM, de- 
pends ultimately on the emission strength relative to the dif- 
fuse background light. This background can be divided into 
four components: terrestrial airglow, zodiacal light, diffuse 
galactic emission (i.e., scattered starlight), and the extragalac- 
tic background. The first is obviously unimportant for space- 
based observations and we will not discuss it further here. The 
other three are unavoidable. 

The best measurements of the backgro unds level in the rel- 
evant regime come from B ernstein et aLl l^OO?). They find to- 
tal backgrounds of 72i ~ 16-50 over the range 3000-8000 A, 
where J^, = J21 x 10~^' erg s~' cm~^ Hz~' sr^'. Of this to- 
tal, > 85% comes from the zodiacal light, which can be mod- 
eled as reflected sunlight. They estimate that diffuse galac- 
tic emission has ^21 ^ 0.5 along clear lines of sight through 
the Galaxy. The extragalactic background (from galaxies 
with y > 23 AB magnitudes) is then J21 = (2.0, 1.3, 1.1) at 
(3000,5500,8000) A, with fractional uncertainties of - 50% 
in each case. These are a few times smaller than the prediction 
of Haardt & Madau (2001) at z = 0. 

Lines with fluxes above these levels would be straightfor- 
ward to detect in a background-limited observation. In our 
units, the integrated background over a velocity width Av is 



^-bkgd = 107: 



21 



Av 



20 km s 



photons s ' cm ^ sr ' . 



(4) 



As we have seen in ^5.41 bright and moderately bright pix- 
els typically have cr ^ 10-50 km s^', at least if we neglect 
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radiative transfer. Thus we expect that features with <i> > 
100 photons s^' cm^^ sr^' should be visible against the zo- 
diacal light, and features with $ > 10 photons s~' cm^^ sr^' 
should be visible against the other backgrounds. This 
emphasizes the importance of high spectral resolution ob- 
servations, because the background increases as the spec- 
tral resolution decreases. The typical line widths are 
0.16(Av/20km s^')[(l +z)/4] A (again neglecting radiative 
transfer effects). We therefore recommend spectral resolution 
be a priority in any observations designed to probe the lowest 
surface brightness objects. 

In fact, so long as the background light is either spectrally 
smooth or can be modeled to high precision, observations can 
in principle extend far below the limits stated above. For 
example, spectral features in the zodiacal light mirror those 
in the solar spectrum and can be modeled fairly well, as 
Bernstein et a n il20()l did. TheLya emission would then ap- 
pear as (weak) fluctuations on the residual background. We 
emphasize that the background need not be spatially uniform, 
so long as the spectru m in each p i xel is either smooth or 
can be modeled. Zaldarriaga et alJ J2004h discuss in some 
detail conceptually similar "foreground-cleaning" techniques 
for low-frequency radio observations (for applications, see 
also lFurlan etto et al. 2004b c). 

6. SOURCE CHARACTERISTICS 
6.1. Luminosity Functions 

We will now consider individual Lya sources. We first di- 
vide the simulation particles into gravitationally bound groups 
using a friends-of-friends group finder with linking length 
equal to 20% of the mean interparticle separation and a min- 
imum group size of 32 particles. The linking is performed 
on the dark matter, with each baryonic particle assigned to its 
nearest dark matter neighbor We compute the Lya luminosity 
of each source as described in ^ 

Figure \^ shows some of the characteristics of groups in 
the Q5 simulation at z = 3. Panel (a) compares the Lya lu- 
minosity from the IGM to that from SF. The diamonds, tri- 
angles, and squares assume self-shielded gas is in CIE, has 
Cq, = 0, and remains optically thin, respectively. In each case 
we have randomly selected 300 well-resolved groups to dis- 
play (i.e., they each have at least 10 times the minimum group 
mass). The solid line shows LLya = ^sfr- For groups with 
no ongoing SF, we set Lsf = 10^^ erg s^' if the group has 
star particles and Lsf = 10^^ erg s^' otherwise.^ The scatter 
increases toward lower luminosities, at least in part because 
accretion and star-formation are more stochastic in halos that 
are not as well-resolved. Clearly, the SFR and the IGM lu- 
minosity are roughly proportional in all treatments, as one 
would expect if the IGM emission comes from gas that will 
soon become available for star formation (provided that the 
gas temperature does not vary strongly between the different 
groups). The wind prescription also forces the ejected mass 
to be proportional to the SFR, so it leads to the same relation. 
Lsf typically exceeds the IGM luminosity except in the most 
optimistic case (which is CIE at z = 3). The emission from 
optically thin gas is typically only 0.01-0.05Lsf; if self- 
shielded gas does not emit, the luminosities of most sources 
would be completely dominated by the emission from the as- 
sociated SF. Note, however, that the symbols on the left in- 
dicate that some groups have sizable Lya luminosities even 

' The latter set is composed of groups near the resolution limit. 
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Fig. 12. — Characteristics of bound groups in the Q5 simulation at z = 3. 
(a): Compares the IGM luminosity to that from stars in each group. If the 
group has a vanishing cun'ent star formation rate (SFR) but contains star 
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hya luminosity (including SF) as a function of group mass, with the solid 
line showing the slope if tLya oc M. The different sets of points correspond 
to our different assumptions about self-shielded gas. In each case we have 
randomly selected 300 well-resolved groups. 



though they are not currently forming stars. Of the groups 
above our mass threshold, 13% have no ongoing SF but do 
have star particles, while another 1.7% have neither SF nor 
star particles. 

FigureEJ? shows the total Lya luminosity as a function of 
group mass. This is nearly independent of our treatment of 
self-shielded gas because SF dominates in most cases, except 
of course for objects without ongoing SF (which are normally 
near the mass threshold). As we approach this threshold, the 
scatter increases rapidly because SF is stochastic in poorly- 
resolved halos. The solid line shows LLya c>c M. This appears 
to be slightly too shallow to match the simulations, but only 
by a small amount. Figure 7 of Springel & Hernquist ( 2003q|) 
shows that, in the mass range we consider here, the SFR is 
also roughly proportional to mass, so this is not surprising. 
For higher mass halos (taken from the D5 simulation, for ex- 
ample), the relation steepens slightly. Interestingly, however, 
the IGM emission declines slightly relative to that from SF 
in massive halos. If we had included only IGM emission, 
Llyq oc M remains a reasonable fit, although the relation is 
slightly steeper if we exclude emission f rom self-shielded gas . 

This contrasts with the predictions of iHaiman et al.l J2000I) 
and Fardal et al. (2001), who argue that cooling following 
gravitational collapse should yield a total luminosity propor- 
tional to the binding energy of the halo, or Ljcm oc M^^^. We 
find that the photoionizing background actually provides most 
of the energy in our calculations; previous work had neglected 
this component. In this case, the Lya luminosity should de- 
pend primarily on the mass of cool gas in each group, so we 
naturally expect LLya oc M. The fraction of energy that can 
be traced to the ionizing background is smallest if = for 
self-shielded gas, which explains why the relation is steep- 
est in that case. Of course, some of this energy exchange is 
unphysical, de pend ing on if and when the gas becomes self- 
shielded (see 44-.21 . If most of the contraction occurs while 
gas is shielded, the only available heat source would be grav- 
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Fig. 13. — The evolution of the luminosity function of groups, (a): Total luminosity function, including both IGM (assuming CIE for self-shielded gas) and 
SF. (b): IGM emission only, assuming CIE for self-shielded gas. (c): IGM emission only, assuming all gas is optically thin, (d): IGM emission only, assuming 
= for self-shielded gas. All distances are comoving. Error bars (shown only for z = 3) are Poissonian. 



itational and we would expect Llyq cx M^/^. The difference 
could also be due to our wind prescription. Winds eject gas 
more efficiently in smaller halos and provide a new source 
(beyond accretion) for the emitting particles. However, al- 
though decreasing the strength of the winds increases Lsf by 
a factor of a few, it has little effect on the IGM emission. 
The IGM luminosity of small groups increases slightly for 
the weak wind case, but it does not appear to be statistically 
significant. 

Figure shows how the comoving luminosity function of 
these groups evolves with redshift. Panel (a) shows the to- 
tal Lya luminosity of each group, while the others include 
only IGM emission. Figures 113b and 1131 ' assume that the 
self-shielded gas is in CIE, Figure [Oj: assumes that all gas 
is optically thin, and Figure [T^ sets = for self-shielded 
gas. We use the Q5 simulation for z = 5, the D5 simulation for 
z = 3 and 2, and the G5 simulation for z = 0. We show sam- 
ple error bars for the D5 simulation at z = 3 assuming Poisson 
statistics. Note that the sharp cutoffs at large luminosities for 
the z > 2 simulations are not real; they are due to the finite 
size of the simulation boxes. (However, the difference across 
panels in the location of the cutoff is real.) The turnover at 
small luminosity for z < 3 is due to the finite mass resolution 
and is also not real. 

The trends with redshift are essentially as expected. As 
emphasized above, the total emission is dominated by SF, so 
panel (a) qualit atively resemb l es the "SF multiplicity func- 
tion" of Springel & Hernauist* ('2003b) (except that we have 
not divided low-redshift groups and clusters into their con- 
stituent galaxies). In particular, dn/dlnL declines at z < 1 



along with the global SFR (Springel & Hernauist'2003b'). The 
IGM luminosity functions depend principally on the treatment 
of self-shielded gas. Groups at z ~ 2^ are most luminous if 
we assume CIE, while those at z and z 5 are most lumi- 
nous if all gas is optically thin. As before, this is a direct con- 
sequence of the temperature of the cooling locus. The max- 
imum IGM luminosities without self-shielded gas are much 
smaller. If this were the most accurate model then it would 
be considerably more difficult to detect the coronae in Lya 
emission. 

Figure[2]shows the luminosity function at z = 3 in several 
different simulations. Here we show the cumulat ive f unction 
for comparison with existing observations (see ^6.3> . Panel 
(a) again shows the total Lya luminosity of each group, while 
panel (b) includes only the IGM emission. In both cases we 
assume that the self-shielded gas is in CIE. These show clearly 
that our luminosity functions have converged over most of the 
appropriate range and that the low and high LLya cutoffs are 
not physical. We find similar behavior for the other treatments 
of dense gas. 

As in 13 we find that the IGM emission depends sensi- 
tively on our assumptions about dense gas: the luminosity 
function of this component varies by several orders of mag- 
nitude in the different cases. Thus observations of (or con- 
straints on) the diffuse emission around galaxies will help to 
determine how gas accretes onto galaxies and how this gas in- 
teracts with the ionizing b ackground (d i ffus e and local) and 
with galactic winds. Like Fardal et al. (2001), we find that 
the Lya emission from SF typically dominates the total lu- 
minosity, even in the most optimistic IGM treatments. How- 
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Fig. 14. — The Lya luminosity function of groups at z = 3 in the different 
simulations, (a): Total luminosity function, including both IGM (assuming 
CIE) and SF. (b): Luminosity function for IGM emission only, assuming GIF. 
All distances are comoving. Error bars are Poissonian. 



ev er, the dispar i ty betw een the two components is not as large 
as iFardal et alJ J2001I) claim unless = for self-shielded 
gas (see their Figure 2). One possible explanation is that the 
conservative formulation of SPH employed in our simulations 
is much less susceptible to numerical overcooling problems 
that trouble d earlier work (Scrinael & Hernauist 2002). An- 
other is that iFardal et"al] 0200 1 ) used a simulation without a 
photoionizing background to compute the IGM emission. We 
have argued that this is not a good assumption in general, al- 
though it may be adequate if the gas becomes self-shielded 
before most of the gravitational heating. In that limit their 
method is in principle a cleaner test of the importance of grav- 
itational shock-heating. Their results lay between our three 
cases, suggesting that (1) a substantial fraction of the gravita- 
tional energy is radiated after the gas becomes self-shielded 
and/or (2) more of the gravitational energy is radiated through 
other cooling mechanisms in our simulations. The former is 
not inconsistent with our results, because the self-shielded gas 
has contracted more than the optically thin component, which 
lies in the outskirts of the coronae. We also note that SF emis- 
sion is much more likely to suffer from dust extinction, which 
would make the IGM emission easier to see along those lines 
of sight. 

6.2. Characteristic Sizes 

We use SExtractor (version 2.3.2. iBertin & Arnout3ll996h 
to identify sources within our maps. We choose the default 
parameter values with the following exceptions: (1) we set 
the deblending parameter to 0.05 to match the choice of M04, 
(2) we do not convolve the image with a filter before identify- 
ing sources, and (3) we identify sources using an absolute sur- 
face brightness rather than a multiple of the local background. 
We make this choice because our maps are (by construction) 
background-free and we wish to avoid any artifacts from im- 
posing one. 

To give some intuition, the green circles in Fig- 
ure |6l outline sources identified by SExtractor with a 
threshold of $ = 100 photons s^' cm~^ sr~^ (or 9.6 x 
10~^' erg s^' cm^^ arcsec^^). In each case, the area of 
the circle is proportional to the isophotal area of the object 
(though we do not outline the shape of the individual sources). 
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Fig. 15. — Scatter plot of angular size against mean surface bright- 
ness for sources from the Q5 simulation at z = 3. Panels (a) and (h) 
have * = 7000 and 100 photons s"' cm-^ sr"' (or 6.7 X 10" and 9.6 X 
10~^' erg s^' cm^- arcsec^-), respectively. The triangles are for emission 
from SF; other symbols are for IGM emission with the different treatments 
of self-shielded gas. Note that 1" = 21.6/i~' comoving kpc. 



We have set the minimum area to 2 square arcseconds for 
completeness here. The largest source in the left slice has 
A = 237 square arcseconds for SF. It splits into two sources 
for the IGM, with A = 185 and 182 square arcseconds. Only 
five (or three) other sources have A > 10 square arcseconds 
for IGM (or SF). In the right-hand map, only two sources have 
A > 10 square arcseconds, in both IGM and SF. In general, the 
IGM Lya emission is more spatially extended but fainter than 
the associated emission from SF. 

Figure ^] shows a scatter plot of isophotal area versus 
mean surface brightness for slices from the Q5 simulation 
at z = 3. We have used ten slices of thickness Az = 0.01 
and width 5.5/;^' Mpc (comoving) to construct the sample. 
The angular resolution is 1" (equal to 21.6/i^' comoving 
kpc). The total volume is thus about twice the volume 
of the Q5 box (but most repeated sources will be viewed 
from different directions). In all cases we have imposed 
a minimum source area of 5 square arcseconds. The tri- 
angles are for SF, while the other symbols show the IGM 
emission for our three cases. Figures I15fe and b have de- 
tection thresholds of $ = 7000 photons s^' cm^- sr^' 
(or 6.7 X lO^'** erg s^' cm~^ arcsec^-) and 
= 100 photons s^' cm^^ sr~' 
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$ = lOOphotonss^ cm^ sr~ (or 9.6 

10^^' erg s^' cm^- arcsec"^), respectively. For the for- 
mer, there are no IGM emitters that meet this threshold 
unless the self-shielded gas is in CIE; the other samples 
are complete. For the lower threshold of Figure 1151 '. we 
randomly select 250 objects from the samples (unless eo, = 
for self-shielded gas, in which case only 16 sources meet our 
criteria). 

Comparing the IGM and SF points at the two $ thresholds, 
we see that the IGM emission is usually fainter but more ex- 
tended than the SF component. This remains true until the 
highest thresholds, when the IGM emission cuts off before SF 
emission; of course the location of the cutoff depends on the 
treatment of dense gas. At the flux levels currently accessible 
to observations (FigureFTSfc) IGM emission will be difficult to 
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Fig. 16. — The distribution of Lyo emitters as a function of isophotal area 
at z = 3. The angular resolution is 1" (21.6/)"' comoving kpc). (a): The 
solid curves show total emission, while the dashed and dot-dashed curves 
show emission from the IGM only, assuming CIE for self-shielded gas and 
that all gas is optically thin, respectively. Thick and thin curves set the thresh- 
old at $ = lO' and 2.3 X lO" photons s~' cm"- sr"' (or 9.6 X 10"-" and 
2.2 X 10~'* erg s~' cm~^ arcsec"^). All curves are from the Q5 simula- 
tion, (b): The solid, dashed, and dotted curves are for total (IGM and star 
formation) emission in the Q5, D5, and G6 simulations. The cutoff at small 
area is due to the SExtractor threshold that we chose. 
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Fig. 17. — The distribution of Lyo emitters as a function of mean surface 
brightness at z = 3. The angular resolution is 1" (21.6/;^' comoving kpc). 
(a): The solid, short-dashed, and long-dashed lines are for the total emission, 
IGM assuming CIE for self-shielded gas, and IGM assuming that all gas is 
optically thin; each of these has A„^i„ = 5 square arcseconds. The dot-dashed 
line shows the total emission with A,„i„ = 16 square arcseconds. (b): The 
solid, long-dashed, short-dashed, and dotted curves are for total emission in 
the Q5, Q3, D5, and G6 simulations, with A,^i„ = 5 square arcseconds. 



detect except where it is self-shielded, has a high temperature, 
and is in CIE. Interestingly, in this case there are about twice 
as many IGM sources with $ > 7000 photons s^' cm^^ sr^' 
as there are SF sources; this is because most of the star- 
forming regions have A < 5 square arcseconds. 

Figures [^and[n]show the source characteristics in a more 
quantitative form. Figure^]shows the angular density of Lya 
sources per unit redshift as a function of isophotal area (note 
that the cutoff at small area is a result of the minimum de- 
tection area we set in SExtractor). Panel (a) compares the 
total (including SF and assuming CIE for self-shielded gas, 
solid curve) and IGM (dashed curves for CIE, dot-dashed 
curve for optically thin gas) emission at two different flux 
thresholds, $ — 10^ and 2.3 x 10"* photons s~' cm^- sr^' (or 
9.6 X 10"^" and 2.2 x 10""* erg s"' cm"^ arcsec"^; thick 
and thin curves, respectively). Note that there are no objects 
with $ > 2.3 X 10"* photons s"' cm^^ sr^' if all gas is opti- 
cally thin. Again, we find that the isophotal area of the total 
emission is dominated by the IGM at small $, but at suffi- 
ciently large surface brightness thresholds, the SF component 
dominat es by a large factor. 

Figure [T6t' shows the total emission for the Q5 (solid), D5 
(dashed), and G6 (dotted) simulations. We also show error 
bars on the Q5 curve, assuming Poisson statistics. Clearly 
the G6 simulation underestimates both the number of faint, 
compact sources and the number of highly extended sources. 
The former occurs because the simulations have a finite mass 
resolution. The latter probably occurs because the larger- 
scale simulations do not fully resolve the spatial distribution 
of SF within halos, even though the global SFR converges 
(see ^5.2> . The higher-resolution simulations do seem to have 
converged reasonably well. 

Figure shows similar information, except as a function 
of the mean surface brightness. Panel (a) shows that the SF 
component dominates the total brightness, especially if the 
dense gas is optically thin (or if it cannot emit). Note also that 
a larger area threshold decreases the abundance of both faint 
and extremely bright sources. FigureFTTb shows poor conver- 
gence in the brightness distribution: the number of sources 
with $ > 10'* photons s"' cm^^ sr^' varies by more than an 
order of magnitude between the simulations. This may be sur- 



prising given that the Lya luminosity function converges rea- 
sonably well. It is again a result of failing to resolve the spatial 
distribution of the gas. The high-resolution simulations bet- 
ter describe the extended, cool gas, so the sources have larger 
areas and hence smaller mean $ (even for the same intrinsic 
luminosity). 

We note that winds have little effect on the sizes of the Lya 
emitters. There is weak evidence that the low surface bright- 
ness IGM emission is slightly more extended in the case of 
strong winds, but it is not a significant difference. This is not 
surprising given that the winds will have difficulty escaping 
along the direction of bright, dense particles. 

In this section we have focused on z = 3, because that is 
where the most complete observations are. However, our 
qualitative conclusions hold at other redshifts as well: star- 
forming regions produce compact, bright Lya emission sur- 
rounded by extended, fainter IGM emission. 

6.3. Comparison to Observations 

In the past several years, Lya line selection has become 
an increasingly popular technique to identify high-redshift 
objects. The usual procedure is to compare broadband and 
narrowband images of a field; those objects that have high 
equivalent widths in the narrowband image are likely to be 
emission-line galaxies. In most cases, contamination from 
other lines (particularly [O II] A3727 at low redshifts) is a 
worry, but color selection and spectroscopic follow-up can 
determine the sample purity. To date, most studies have 
focused on relatively high-surface brightness emission from 
star-forming galaxies. 

A number of these surveys have constrained the space 
density of luminous, relatively compact Lya emitters. We 
will begin by comparing these observations to our results. 
ICowie & Hul ( fT9 98) performed blank-field surveys at z w 3.4 
and found « ~ (3, 1) x IQ^^hr' Mpc^"* for objects with L42 > 
(2,5). Here Llyq ~ L42 x 10^^ erg s~'. This is about an or- 
der of magnitude smaller than the predictions of Figure 
(Note that the luminosity function evolves only weakly with 
redshift in the range z ^ 2-5.) SO O found a space density 
about 6 times larger than ICowie & Hu ( 1998 ); however, the 
SOO volume contains ^ 6 times as many continuum-selected 
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galaxies as a n average field, so this is probably due to cos- 
mic variance. lFuiita etalJ(l20()l find « - 2 x lO^V Mpc"^ 
for objects with L42 > 5 at z = 3.7, several times smaller than 
[Cowie & Hu ( 1998). Ouchi et al. (2003) performed a survey 
for Lya emitters at z = 4.86 and found no evidence for evolu- 
tion in the lum inosity function between that time and z = 3.4. 
iPalunaset all (12004) foundn - 3 x 10"^/i^ Mpc^^ forL42 > 4 
at z = 2.38, also consistent with weak redshift evolution. 

All of these surveys have different selection criteria, so de- 
tailed comparisons to our luminosity function are difficult. 
However, it seems clear that our predictions are roughly an 
order of magnitude above the observed number densities at 
L42 ^ 1-10. There are several possible explanations for the 
discrepancy. First, we have included the total emission from 
SF and the IGM. The IGM emission is relatively extended 
and may not actually contribute to the observed high-surface 
brightness emission. If the CIE case is overly optimistic, this 
would also decrease the luminosities by about 50%. Second, 
we have ignored processes within the galaxies that can dra- 
matically modify the Lya line luminosity, such as dust (which 
appears to completely eliminate or substantially reduc e Lya 
emission in ~ 75-80% of z = 3 galaxies; SOO, Shaple v et al.l 
I2003J . Such a correction would bring our estimates into rea- 
sonable agreement with most of the observations described 
above. The dust extinction probably depends on the line of 
sight through the galaxy (and its wind); much of the galaxy 
will therefore show little Lya emission, allowing the IGM to 
stand out more clearly. Th ird, we have used a crude conver- 
sion of SFR to Lhya (see ^3.3> . A more sophisticated treat- 
ment may help to reconcile the predictions with the observed 
number density. 

To probe the IGM emission, it is best to consider ex- 
tended objects with relatively low surface brightness. SOO 
and M04 have detected extended "blobs" of Lya emis- 
sion without associated (significant) UV continua near a 
"proto-cluster" of Lyman-break galaxies. SOO detected two 
blobs with Li,ya ^ lO'^'^ergs"' and physical extents of ^ 
100/;^' kpc, implying n ~ 3 x Mpc^^ for these ob- 

jects. (Note, however, that this is the same field cited 
above and has ^ 6 times as many galaxies as average). 
M04 surveyed an overlapping field (with effective volume 
about eight times larger) and selected objects with A > 16 
square arcseconds and <i> > 2.3 x 10"* photons s~' cm~^ sr~' 
(2.2 X 10^'^ erg s^' cm arcsec"^); they find n ^ 1 x 
IQ^'^h^ Mpc^^. These objects have L42 > 6. 

Extended emission is thus somewhat rarer than the compact 
Lya emitters probed by most surveys, and our models can 
easily accommodate the raw luminosities of the M04 sample. 
Figure[T4t' shows that the IGM component on its own can pro- 
vide enough energy so long as we assume CIE. However, most 
of the sources in the simulation are much smaller than the 
M04 objects. With their area and surface brightness thresh- 
olds, we find - (3,20) x IQ-'^h'^ Mpc"^ for IGM and 
total emission, respectively, assuming CIE for self-shielded 
gas. In the other cases no IGM objects satisfy these criteria. 
Unfortunately, these predictions are based on only a few ob- 
jects and have not fully converged for the IGM emission. The 
strongest statement we can make is that star formation appears 
to produce easily enough Lya emission on moderately large 
(A ~ 20-40 square arcsecond) scales to explain the observa- 
tions, even if we renormalize the number densities to match 
that of Lya galaxies (thus crudely accounting for dust extinc- 
tion in the SF component). However, the IGM on its own 



could only account for the blobs in the most optimistic case. 

If we consider the most luminous and most extended ob- 
jects, the IGM looks to be even less plausible. At the M04 
$ threshold, none of our objects have A > 60 square arcsec- 
onds even if we include star formation, while ^ 10% of the 
observed blobs exceed this size (and the largest has A — 222 
square arcseconds). We most likely therefore need to appeal 
to other processes to explain the largest and most luminous 
blobs (including those first discovered by SOO), although we 
cannot rule out that the finite size of our simulation box is 
responsible for the difference. We must also stress that the 
observed luminous Lya bl obs have velocity w idths of sev- 
eral hundred km s^' (SOO; iBower et al.ll2004l), m uch larger 
than the intrinsic widths of our sources (see ^5.4> . If these 
sources are ultimately due to cooling gas and SF, large in- 
ternal velocity gradients (from winds, for example) or com- 
plex radiative transfer (e.g., Zheng & Miralda-Escude 200^ 
would have to dramatically alter our predicted linewidths. On 
the other hand, both SOO and M04 selected a region known to 
contain a strong galaxy overdensity, and the density of blobs 
in more normal environments is unknown. M04 note that the 
blobs are highly clustered even within their region. Observa- 
tions of more representative regions would be invaluable in 
understanding the nature of strong Lya emission. 

Another interesting issue is the correlation between IGM 
emitters and galaxies. In our simulations, Lya sources with- 
out associated stellar particles are extremely rare. In most 
cases, those that do exist are associated with poorly-resolved 
halos and may be simulation artifacts. Thus, within the con- 
text of our simulations, isolated bright Lya emitters will be 
rare. Although the luminous SOO blobs are not centered on 
optically visible galaxies, at least one of them has a vigor- 
ously star-forming (SFR > 5OO/r_^M0 yr^') submm galaxy 
near its center (jChapman et al.l l2001h and both have optical 
galaxies near their outskirts. M04 found that 7/33 other blobs 
match with previously detected galaxies. The others may also 
be associated with faint galaxies (perhaps due to intrinsic ab- 
sorption), and testing this connection would be useful in the 
future. Such observations could also help to measure the rel- 
ative importance of Lya from SF and the IGM in different 
galaxies. We have found that SF usually dominates the to- 
tal luminosity of each system; however, we have also found 
that our simulations overestimate the number density of Lya 
-emitting galaxies. This is likely due to internal structure 
and the dust distribution within galaxies. The IGM emission 
should be less susceptible to these uncertainties, so in some 
cases it may actually dominate the observed emission. 

Finally, we have neglected radiative transfer in computing 
the source sizes. Because our regions are optically thick in 
the Lya line, these photons scatter resonantly until reaching 
an optical depth near unity. Although most of the scattering 
occurs in frequency space, some photons could in principle 
also scatter in real space to the edge of the optically thick 
region, from which they can escape even from the center of 
the line. Thus our method may underestimate the true ex- 
tent of the emitting regions. Fardal et al. (2001) argue that 
this mechanism can help to reconcile the luminous, compact 
sources found in their simulations with the significantly more 
extended observed systems. Furthermore, ionizing photons 
from young stars in the central galaxy can escape into the sur- 
rounding IGM but get absorbed by the cool, dense gas. This 
would shift some of the Lya emission that we attribute to star 
formation into the IGM, making the SF sources larger than we 
have found. 
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7. DISCUSSION 

We have used cosmological simulations to study the Lya 
emission expected from structure formation over a broad 
range of redshifts. We considered emission from optically 
thin low-density gas, from self-shielded gas clouds, and from 
the reprocessing of stellar ionizing photons. As in F03, who 
studied Lya emission from z < 0.5, we have found that the 
brightest regions are compact star-forming galaxies. These 
are surrounded by extended coronae of emission from IGM 
gas; the regions closest to galaxies are likely to be self- 
shielded (at least from the mean ionizing background) while 
the outer reaches remain optically thin. Crucially, although 
emission associated with star formation dominates the total 
Lya luminosity, the IGM coronae are more extended and can 
thus be isolated (at least in principle). They could also dom- 
inate in galaxies in which dust has extinguished the Lya line 
from star formation. However, the luminosities of the coronae 
depend sensitively on our assumptions about gas that is able 
to shield itself from the metagalactic ionizing background. 

One of the most important results of this study is that the 
photoionizing background can power the bright Lya coronae, 
in contrast t o earlier predictions i n which gravitatio nal energy 
dominated (Haiman et al. 2000; Fardal et al.l 1200 1). In our 
simulations, most of the flux comes from gas on a high den- 
sity (hh ^ lO^-' cm^-'), low temperature (T < lO"^-^ K) track 
of particles bound to dark matter halos but outside of galaxies. 
Although many of these dense particles have been moderately 
shocked in the past (see Figure |5}, the post-shock cooling 
times are so short that they quickly enter a quasi-equilibrium 
cool phase. In this regime, photoheating and radiative cooling 
nearly balance, increasing the cooling time. While it is opti- 
cally thin, the ionizing background powers > 50% of the Lya 
emission at z > 1, provided that it is near or above our fiducial 
amplitude. As the gas condenses, we expect it to eventually 
become self-shielded from the diffuse ionizing background. 
Because we cannot include radiative transfer in our calcula- 
tions, the power source beyond this point is unclear. If grav- 
itational contraction or winds maintain the gas temperature 
near the peak for Lya excitation, photoionization would play 
a relatively minor role. If, on the other hand, the dense gas 
rapidly cools or remains optically thin (in particular due to 
embedded sources), photoionization can power > 90% of the 
Lya emissivity. We thus argue that including photoionizing 
radiation is crucial to interpreting observations of Lya cool- 
ing around galaxies. 

In any of these cases, the cool gas in Lya coronae directly 
traces the accretion of gas onto galaxies, although winds do 
significantly complicate the interpretation. In the presence of 
winds, Lya emission comes from a mix of relatively pristine 
accreting gas and material that has recently been ejected from 
the central galaxy but remains trapped inside the halo (either 
because the wind velocity is smaller than the escape speed 
of the halo, or because the wind interacts with the accreting 
gas). We find, however, that winds have a surprisingly small 
effect on the Lya emission (especially considering that they 
reduce the global SFR by a factor of five), although our sim- 
ulations do not fully resolve the winds. Clumping within the 
wind medium could increase the local emissivity, or the faster 
cooling may end up reducing the total emission. 

Thus the Lya coronae present the intriguing opportunity to 
measure some combination of wind feedback near and gas 
accretion onto galaxies. Unfortunately, the self-shielded gas 
responsible for most of the emission is also the most diffi- 



cult to model. We have emphasized several uncertainties in 
^3.21 including the presence of dust, the geometry and veloc- 
ity structure of the region, radiative transfer of higher Lyman 
line photons through the cloud, the temperature of the self- 
shielded gas, and the local ionizing radiation field. Of these, 
the last two are likely the most important. 

The temperature of gas on the cooling locus is determined 
in the simulation by the balance between photoheating and ra- 
diative cooling. If the cloud is self-shielded, the primary heat 
source disappears and it may cool rapidly. The Lya emis- 
sivity is extremely sensitive to the temperature in collision- 
ally ionized gas, so small errors in the heat balance could 
translate into large changes in the emission properties. Un- 
fortunately, the simulations do not include self-shielding be- 
cause they do not incorporate a self-consistent treatment of 
radiative transfer Although we can crudely identify regions 
that are shielded from the metagalactic ionizing background, 
we cannot account for radiation from nearby galaxies or self- 
consistently describe the behavior of shielded gas. We there- 
fore consider three cases that should bracket the real behavior 
In the first, we set the emissivity of self-shielded gas to zero, 
which would be a good approximation if the gas were to cool 
efficiently. In the second, we assume that the simulation tem- 
peratures are accurate and that the self-shielded gas is in CIE. 
In the third, we assume that the local ionizing field is much 
larger than the mean, so that even dense gas remains opti- 
cally thin. In the first case, IGM gas is extremely faint and 
may be beyond the reach of near-future surveys. But in ei- 
ther of the other cases, the gas is significantly brighter and the 
total luminosity is not too far below that from star formation. 
The balance between the two optimistic cases depends on red- 
shift, because the emissivity in CIE is extremely sensitive to 
temperature, which, in our simulations, is set by the ionizing 
background. 

We have also computed the Lya emission from recombina- 
tions following absorption of stellar ionizing photons inside 
galaxies. We have used an extremely simp le co nversion be- 
tween local SFR and Lya luminosity (see ^3.3> . appropriate 
for a standard IMF and little extinction. The true luminos- 
ity will vary with the initial mass function, metallicity, dust 
content, and internal structure of each galaxy, so our quan- 
titative results should be viewed as no more than represen- 
tative. With these assumptions, this component dominates 
both the total luminosity and the brightest pixels, but it is 
confined to smaller physical scales than the IGM emission. 
We have, however, assumed that all of the ionizing photons 
are absorbed locally. In reality, some fraction of the ionizing 
photon s will escape the host galaxy into the surrounding IGM 
( Steid el et al.ll200U Eeaueux et al. 1995; Kunth et al. 2003|. 
It is this radiation that could escape and help to illuminate the 
nominally self-shielded gas in the Lya coronae. By neglect- 
ing dust, we have also overestimated the total amount of SF 
Lya radiation by a factor of ~ 4 (at z = 3; SOO). Both of 
these effects help to increase the relative importance of IGM 
emission. 

A number of recent observations have targeted high- 
redshift galaxies with strong Lya emission (see ^6.3> . The 
space density of such objects is about an order of magni- 
tude smaller than predicted by our simulations. A large part 
of this difference is likely because only ~ 20-25% of star- 
forming galaxies have strong Lya emission (SOO), for the rea- 
sons outlined above. The remaining discrepancy is reason- 
able given the simple assumptions behind our model. Our 
results therefore suggest that Lya selection can be a useful 
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probe of galaxy formation and gas accretion, even if high- 
surface brightness galaxies are the only targets. This is par- 
ticularly interesting given that Lya -selection is currently the 
most e fficient technique to select the highest-redshift galaxies 
(e.g., Hu et al. 2002; Kodaira et al. 2003; Rhoads et al. 2004 
and can pot entially be a powerful probe at even earlier epochs 
(^arton et al.l2004l; though absorption from the neutral inter- 
vening IGM could complicate such an interpretation; see, e.g., 
iFurlanetto et al. 2004d and references therein). 

We have emphasized the substantial uncertainties in Lya 
emission from both the IGM and star formation. While 
these imply that our predictions are far from rigorous, 
they present an excellent opportunity to constrain the 
models through observations. Existing surveys for diffuse 
Lya emission at z 2-3 have reached surface bright- 
ness thresholds of ~ 2.3 x 10"* photons s^' cm^^ sr^' (or 
2.2 X 10""^ erg s"' cm'^ arcsec"^ at z = 3.1; M04) with 
the Subaru telescope and ~ 1.3 x lO'* photons s^' cm^^ sr^' 
(or 1.6 X 10~'^ erg s~' cm~^ arcsec~^ at z = 2.2; 
[Francis & Bland-Hawthorn' •2007) with the 4-meter Anglo- 
Australian Telescope. Both are fairly close to our more 
optimistic predictions. We have shown that the combi- 
nation of star formation and cooling gas powered by the 
photoionizing background can in principle account for most 
of the Lya blobs observed by M04, although we found no 
counterparts to the largest and most luminous blobs identified 
by SOO. IGM emission alone probably cannot explain the 
observations, but in our more optimistic cases should appear 
with order-of-magnitude increases in the surface brightness 
sensitivity. 

We find no evidence in our simulations for strong Lya 
emission except that surrounding galaxies. We therefore ex- 
pect that all the Lya blobs observed at z 3 contain em- 
bedded galaxies, which could be tested through follow-up 
observations of the M04 sources. An alternate approach 
to find isolated emitters is t o seek quas a r abso rption line 
systems in emission (iHogan & Wevmannlll987h . Analytic 



estimates suggest that optically-thick Lyman-limit systems 
should be visible through reproce ssing of the ionizing back- 
ground (iGould & Weinl3erg|[T996l) . In our simulations, it is 
most natural to associate these with the Lya coronae around 
galaxies, so it is important to identify the absorber envi- 
ronments through observations. Such searches are already 
underway. Francis & Bland-Hawthorn (2004) searched for 
diffuse emission near a bright quasar at z = 2.168. They 
found no sources with A > 9 square arcseconds and $ > 
L3 X 10"* photons s^' cm^^ sr^'. This is somewhat surpris- 
ing in our model, given that we expect a quasar to be embed- 
ded inside an overdensity with an abundance of cool gas; the 
peak brightness should be at least comparable to their thresh- 
old if the gas emits in collisional ionization equilibrium. The 
most likely explanation is that the quasar has photoionized 
all of the nominally self-shielded gas. This would reduce the 
peak to $ ~ 10^-10"^ photons s~' cm~^ sr~' (see FigurelSj'), 
just below the existing threshold. We therefore argue that it is 
crucial to continue to push these na r rowb and observations for- 
ward. Francis & Bland-Hawthorn (2004) have demonstrated 
that valuable upper limits can be obtained even from 4-meter 
class telescopes, and our results suggest that direct detection 
of diffuse Lya emission may not be far away. 
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